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ABSTRACT 



Aims. Photodissociation by UV light is an important destruction mechanism for carbon monoxide (CO) in many astrophysical envi- 
ronments, ranging from interstellar clouds to protoplanetary disks. The aim of this work is to gain a better understanding of the depth 
dependence and isotope-selective nature of this process. 

Methods. We present a photodissociation model based on recent spectroscopic data from the literature, which allows us to com- 
pute depth-dependent and isotope-selective photodissociation rates at higher accuracy than in previous work. The model includes 
self-shielding, mutual shielding and shielding by atomic and molecular hydrogen, and it is the first such model to include the rare 
isotopologues C l7 and 13 C 17 0. We couple it to a simple chemical network to analyse CO abundances in diffuse and translucent 
clouds, photon-dominated regions, and circumstellar disks. 

Results. The photodissociation rate in the unattenuated interstellar radiation field is 2.6 x 10" 10 s , 30% higher than currently adopted 
values. Increasing the excitation temperature or the Doppler width can reduce the photodissociation rates and the isotopic selectivity 
by as much as a factor of three for temperatures above 100 K. The model reproduces column densities observed towards diffuse clouds 
and PDRs, and it offers an explanation for both the enhanced and the reduced jV( 12 CO)/iV( CO) ratios seen in diffuse clouds. The 
photodissociation of C 17 and 13 C 17 shows almost exactly the same depth dependence as that of C 18 and 13 C 18 0, respectively, so 
17 and 18 are equally fractionated with respect to 16 0. This supports the recent hypothesis that CO photodissociation in the solar 
nebula is responsible for the anomalous 17 and 18 abundances in meteorites. Grain growth in circumstellar disks can enhance the 
/V( 12 CO)//V(C 17 0) and 7V( 12 CO)/A r (C 18 0) ratios by a factor of ten relative to the initial isotopic abundances. 

Key words, astrochemistry - molecular processes - molecular data - ISM: molecules - ISM: clouds - planetary systems: protoplan- 
etary disks 



1. Introduction 

Carbon monoxide (CO) is one of the most important mole- 
cules in astronomy. It is second in abundance only to molec- 
ular hydrogen (H2) and it is the main gas-phase reservoir of 
interstellar carbon. Because it is readily detectable and chem- 
ically stable, CO and its less abundant isotopologues are the 
main tracers of the gas properties, structure and kinematics in 
a wide vari ety of astrophysical environments (for recent ex- 
amples, see iDame et al.ll200U iNaiita et al.ll2003t IWilson et afj 



Brown & Miliar! | 2003b iT evssier et all 12004 ICernicharol [20041: 
Morata & Herbstl 120 08). CO ice on the surfaces of grains can 



2005; Greve et al. 2005; Leroy et al. 2005; Huggins et al. 2005; 
Bave t et al l 120061; lOka et al l 120071; [Narayanan et all l2008l) . In 



particular, the pure rotational lines at millimetre wavelengths are 
often used to determine the total gas mass. This requires knowl- 
edge of the CO-H2 abundance ratio, which may dif fer by sev- 
eral orders of magnitude from one object t o the next (lLacv et al.l 
ll994tlBurgh et alj2007tFanic et al.ll2008l) . If isotopologue lines 
are used, the isotopic ratio enters as an additional unknown. 

CO also controls much of the chemistry in the gas phase 
and on grain surfaces, and is a precursor to more complex mol- 
ecules. In photon-dominated regions (PDRs), dark cores and 
shells around evolved stars, the amount of carbon locked up 
in CO compared with that in atomic C and C + determines the 
abundances of small and large carbon-chain molecules (e.g. , 
MillaretaDEHl Uansen et al.lfl99l lAikawa & Herbstl fl999: 



be hydrogenat ed to more complex saturated molecules such a s 
CH3OH (e.g.. lCharnlev et alJll995t[Watanabe & Kouchill2002l) . 
so the partioning of CO between the gas and grains is im- 
portant for the overall chemical composition as well (e.g. , 
Casel li et al. ll!993HRodgers & Charnlevll2003r. iDotv etai]l2004t 
barrod & Herbstll2006l) " 

A key process in controlling the gas-phase abundance 
of 12 CO and its isotopologues is photodissociation by ul- 
traviolet (UV) photons. This is governed entirely by dis- 
crete absorptions into predissociative excited states; any pos- 
sible contributions from continuum c hannels are negligi- 
ble (iHudsonl [l97Tt iFocketalJ [19801 iLetzelter et all fl987[ 
Cooper & Kirbvl Il987l) . Spectroscopic measurements in the 
laboratory at increasingly higher spectral resolution have 
made it possib le for detailed photodissoci a tion models to 
be co n structed (|Solomon &" Klempere rl 1 19721: iBallv & L anger 



1982; Glassgold et al. 1985; van Dishoeck & Black 198f 



Viala et al. 1988; van Dishoeck & Black 1988, hereafter vDB88; 
Warin et all 19961;" Lee et al. 19961) . The currently adopted pho- 



todissociation rate in the unattenuated interstellar radiation field 
is2x 10- 10 s _1 . 

Because the photodissociation of CO is a line process, it is 
subject to self-shielding: the lines become saturated at a I2 CO 
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column depth of abo ut 10 15 cm' 2 , and the pho to dissociation rate 
strong ly decreases dvDB88t iLee et al.l 119961) . iBallv & Langeri 
(1982) realised this is an isotope-selective effect. Due to their 
lower abundance, isotopologues other than 12 CO are not self- 
shielded until much deeper into a cloud or other object. This 
results in a zone where the abundances of these isotopologues 
are reduced with respect to 12 CO, and the abundances of atomic 
13 C, 17 and 18 are enhanced with respect to 12 C and le O. 
For example, the C n O- 12 CO and C 18 0- 12 CO column density 
ratios towards X Per are a fa ctor of five lower than the ele- 
mental oxygen isotope ratios dSheffer et alj|2002l) . The 13 CO- 
12 CO ratio along the same line of sight is unchanged from the 
elemental carbon isotope ratio, indicating that 13 CO is replen- 
ished through low-temperature isotope-exchange reactions. A 
much larger sample of sources shows jV( 13 CO)//V( 12 CO) column 
density ratios both enhanced and reduced by up to a factor of 
two relative to the elemental isotopic ratio dSonnentrucker et all 
l2007tlBurgh et alJfcOOTlfSheffer et alj|2007l) . The reduced ratios 
have so far defied explanation, as all models predict that isotope- 
exchange reactions prevail over selective photodissociation in 
translucent clouds. 

CO self-shielding has been suggested as an explanation 
for th e anomalous 17 Q- 18 abundance ratio found in mete- 
orites dClavton et alJl97Hciavtonll2002l:lLvons & Youngl200l 
ILee et al. 1 120081) . In cold environments, molecules such as wa- 
ter (H2O) may be enhanced in heavy isotopes. This so-called 
isotope fractionation process is due to the difference in vibra- 
tional energies of Hi e O, H^O and H^O, and is therefore mass- 
dependent. It results in ls O being about twice as fractionated 
as n O. However, n O and 18 are nearly equally fractionated 
in the most refractory phases in meteorites (calcium-aluminium- 
rich inclusions, or CAIs), hinting at a mass-independent frac- 
tionation mechanism. Isotope-selective photodissociation of CO 
in the surface of the early circumsolar disk is such a mecha- 
nism, because it depends on the relative abundances of the iso- 
topologues and the mutual overlap of absorption lines, rather 
than on the mass of the isotopologues. The enhanced amounts 
of 17 and 18 are subsequently transported to the planet- and 
comet-forming zones and eventually incorporated into CAIs. 
Recent observations of 12 CO, C 17 and C 18 in two young 
stellar objects support the hypothesis of CO photodissociation 
as the cause of th e anomalous oxygen isotope ratios in CAIs 
( Smi th et alj2009t) . A crucial point in the lLvons~& Young model 
is the assumption that the photodissociation rates of C n O and 
C 18 are equal. Our model can test this at least partially. 

Detailed descriptions of the CO photodissociation process 
are also important in other astronomical contexts. The circum- 
stellar envelopes of evolved stars are widely observed through 
CO emission lines. The measurable sizes of these envelopes are 
limited primarily by the photodis sociation of CO in the radia- 
tion field of background starlight (Mamon et al. I ll988h . Finally, 
proper treatment of the line-by-line contributions to the pho- 
todissociation of CO may affect the anal ysis of CO photoche m- 
istry in the upper atmospheres of planets dFox & Blacklll989b . 

In this paper, we present a n updated version of the photodis- 
sociation model from vDB88, based on laboratory experiments 
performed in the past twenty years (Sect. [2}. We expand the 
model to include C n O and 13 C 17 and we cover a broader range 
of CO excitation temperatures and Doppler widths (S ects. [3]and 
|4|. We rederive the shielding functions from vDB88 and extend 
these also to higher excitation temperatures and larger Doppler 
widths (Sect. [5]). Finally, we couple the model to a chemical net- 
work and discuss the implications for translucent clouds, PDRs 



and circumstellar disks, with a special focus on the meteoritic 
ls O anomaly (Sect.[6]i. 



2. Molecular data 

The photodissociation of CO by interstellar radiation occurs 
through discrete absorp tions i n to pre dissociated bound states, 
as first su ggested by iHudsonl d 197 lb and later confirmed by 
iFock et aiTdl980l) . Any possible contributions from continuum 
channels are negligible a t wavelengths longer than the Lyman 
limit of atomic hydrogen dLetzelter et al.lfl987l : ICooper & K irbv 
fl987h . 

Ground-state CO has a dissociation energy of 1 1.09 eV and 
the general interstellar radiation field is cut off at 13.6 eV, so 
knowledge of all absorption lines within that range (911.75- 
1117.80 A) is required to compute the photodissociation rate. 
These data were only partially available in 1988, but ongoing 
laboratory work has filled in a lot of gaps. Measurements have 
also been extended to include CO isotopologues, providing more 
accurate values than can be obtained from theoretical isotopic 
relations. Table[TJlists the values we adopt for 12 CO. 



2. 1. Band positions and identifications 



Eidelsber g"& Rostasl (I1990L hereafter ER90) and 
Eidelsber g et al. ( 1992I) redid the experiments of iLetzelter et al.l 
d!987l) at higher spectral resolution and higher accuracy, and 
also for 13 CO, C ls O and 13 C I8 0. They reported 46 predisso- 
ciative absorption bands between 11.09 and 13.6 eV, many of 
which were rotationally resolved. Nine of these have a cross 
section too low to contribute significantly to the overall dissoci- 
ation rate. The rem aining 37 bands are largely the same as the 
33 bands of vDB88; bands 1 and 2 of the latter are resolved into 
four and two individual bands, respectively. Throughout this 
work, band numbers refer to our numbering scheme (Table [T), 
unless noted otherwise. 

Th anks to the higher resolution and the isotopologue data, 
ER90 could identify the electronic and vibrationa l charact er of 
the upper states more reliably than Stark et al. in vDB88l The 
vibrational levels are required to compute the positions for those 
isotopologue bands that have not been measured directly. Nine 
of the vDB88 bands (not counting the previously unresolved 
bands 1 and 2) have a revised v' value. 

The lEidelsberg et all d 19921) positions (vo or Ao) are the best 
available for most bands, with an estimated accuracy of 0.1- 
0.5 cm -1 . Seven of their 12 CO bands were too weak or dif- 
fuse for a reliable analysis, so their positions ar e accurate only 
to within 5 cm -1 . Nevertheless, we adopt the lEidelsberg et akl 
positions for three of these: bands 2A , 6 and 14. The former 
was blended with band 2B in vDB88, and the other two show 
a bett er match with the isotopologue band positions if we take 
the lEidelsberg et al. values. For the other fo ur weak or diffuse 
bands, Nos. 4, 15, 19 and 28, we keep the vDB88 positions. 
lUbachs et alj (1994) further improved the experiments, obtain- 
ing an accuracy of about 0.01-0. 1 cm -1 , so we adopt their band 
positions where available. Finally, we adopt the even more accu- 
rate positions ( 0.003 cm -1 or better) available for the CI, £0, El 
and LP bands dUbachs eTaLI l200d: ICacciani et all 1200 ll 120021: 
Cacciani & Ubachs 2004)Q 



1 All transitions in our model arise from the v"=0 level of the elec- 
tronic ground state. We use a shorthand that only identifies the upper 
state, with CI indicating the C 'S + v' = l state, etc. 
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Band positions for isotopologues other than CO are still 
scarce, although many more are currently known from exper- 
iments than in 1988. The CI and £1 bands have been mea- 
sured for all six natural isotopologues, an d the £0 band for all 
but 13 C 17 Q, a t an ac c uracy of 0.003 cm"' (ICacciani et al.ll 19951; 
Ubachs et alj 120001; ICacciani et ail 120011; ICacciani & Ubachsl 
20041) . The positions of the £0 band are especially important 
because of its key role in the isotope-selective nature of the CO 
photodissociation (Sect. 13.3b . Positions are known at lower ac- 
curacy (0.003-0.5 c m-') for an additional 25 C 18 Q, 30 13 CO 
and 9 13 C 18 Q bands (lEidelsberg etalJI 1992b lUbachs etal.lll994t 
ICacciani et al I I2OO2I) : these are included throughout. 

We compute the remaining band positions from theoretical 
isotopic relations. For band b of isotopologue i, the position is 

v (b, i) = v (£, 12 CO) + [E' v (b, i) - E' V '(X, 0] - 

[E' v (b, 12 CO) - E' y '(X, 12 CO)] , (1) 

with E' w (b, i) and E' V '(X, i) the vibrational energy of the excited 
and ground states, respectively. Hence, we need the vibrational 
constants (aj' e , a> e x' e , . . . ) for all the excited states other than the 
C Thes e have only been determined experimentally for the 
E lr l state (Kepa 1988). For the other states, we employ this 
scheme: 

• if it is part of a vibrational series (such as band 30, for which 
the corresponding v'-l band is No. 27), we can derive u)' e 
from the difference in vo; 

• else, if it is part of a Rydberg series convergi ng to the X 2 £ + 
or A 2 n state of CO + , we take those constants (Haridas s et alj 
120001) : 

• else, we take the constants of ground-state CO 
dGuelachvirietal]|1983l) . 

The choice for each band and the values of the constants are 
given in Table Q] 

2.2. Rotational constants 

The rotational constants (B' v and D' v ) for each excited state 
are needed to compute the positions of the individual absorp- 
tion lines. ER90 provided B' v values for most bands, at an es- 
timated accuracy of better than 1%. Their D' v values are less 
well constrained and may be off by more than a factor of 
two. However, this is of little importance for the low-/ lines 
typically involved in the photodissociation of CO. More ac- 
curate values (B' v to better than 0.1%, D' v to 10% or better) 
are available for 12 states from higher-resolution experiments 
dEikema et al] 11994 lUbachs et alj 1 1994 l2000t ICacciani et ail 
I2001L I2002t ICacciani & UbachslfcOOlh Again, the data for iso- 
topologues other than 12 CO are generally scarce, so we have to 
compute their constants from theoretical isotopic relations. This 
increases the uncertainty in B' v to a few per cent. In case of bands 
12, 15, 19 and 28,[ER90 reported constants for C ls O but not for 
12 CO, so we emp loy the theoretical relations for the latter. We 
adopt the vDB88 constants for bands 4 and 14, because they are 
more accurate than those of ER9(J No constants are available 
for Rydberg bands 2A and 6, so we adopt the constants of the 
associated CO + states (A 2 Tl and X 2 1, + , respectively). 

In seven cases, the rotational constants of the P and R 
branch (e parity ) were found to differ from those of the Q 
branch (/parit y; Ubachs et alJl 1994 120001; ICacciani et al1l2002E 
ICacciani & Ub achs 2004). For these bands, the / parity values 
are given in Table Q] Table [2] lists the difference between the / 
and e values, defined as q' v = B' v e — B' y , and p' v = D' ve - D' y ,. 



Table 2. Parity-dependent rotational constants for 12 CO. 



Band 


9t 


P? 


Refs." 


# 


(cnr 1 ) 


(cm 4 ) 




16 


-2.7(-3) 




1 


22 


2.212(-2) 


7.9(-6) 


2 


23 


3.0(-2) 




1 


25 


-7.0(-4) 




1 


26 


-1.1K-2) 




1 


31 


1.14(-2) 


3.0(-8) 


3 


33 


1.196(-2) 


2.1(-7) 


4 



a Dashes indicate that no measurement is available, so we adopt a 
valu e of zero. 

b (1) lUbachs etalJ dl994h: (2) ICacciani et all i2002h : (3) 
lUbachs et all il2000); (4) Caccia ni & Ubachsl ( 120040 



The uncertainties in q' w and p' v are on the order of 1 and 10%, 
respectively. 

2.3. Oscillator strengths 

ER90 measured the integrated absorption cross sections (cri nt ) 
for all their bands to a typical accuracy of 20%, but they cau- 
tioned that some values, especially for mutually overlapping 
bands, may be off by up to a factor of two. The o scillator 
strengths (/ v <o) derived from these data are different from vDB88 
for most bands, sometimes by even more than a factor of two. 
In addition, there are differences of up to an order of magni- 
tude between the cross sections of I2 CO and those of the other 
isotopologues for many bands shortwards of 990 A. The iso- 
topic differences are likely due in part to the difficulty in deter- 
mining individual cross sections for strongly overlapping bands, 
but isotope-selective oscillator strengths in general are not unex- 
pected. For example, they were also observed recently in high- 
resolution measurements of N2 (G. Stark, priv. comm.). For CO, 
the oscillator strengths depend on the details of the interactions 
between the J', v' levels of the excited states and other rovi- 
bronic levels. These interactions, in turn, depend on the energy 
levels, which are different between the isotopologues. We adopt 
the isotope-selective oscillator strengths where available. In case 
of transitions where no isotopic data exist, we choose to take the 
value of the isotopologue nearest in mass. This gives the clos- 
est match in energy levels and should, in general, also give the 
closest match in oscillator strengths, which to first order are de- 
termined by the Franck-Condon factors. 

For ten of our bands, we adopt oscillator strenghts from 
studies that aimed specifically at measuring this parameter 
(lFedermanetal.ll200l[ lEidelsberg et"ai]|2004 50061) . Their es- 
imated acc uracy is 5-15% . The o scillator strength for the £0 
band from Fe derman et alj (|2001) is almost twice as large as 
that of vDB88 and IER901 which appears to be due to an in- 
adequate treatment of saturation effects in the older work. The 
2001 value corr esponds well to other values derived since 1990. 
iFederman et al.l also measured the oscillator strength of the 
weaker CI band and fo und it to be the sa me, within the er- 
ror margins, as that of lvDB88l and lER90t Recent measure- 
ments of the KQ, LQ and £1 transitions and the four W transi- 
tions show larger oscillator streng ths than those of vDB88| and 
ER90 (Eidelsb ergetai]2 004. 2006j). The new values correspond 
closely to those of She ffer et al.1 (120031) . who derived oscillator 
strengths for eight bands by fitting a synthetic spectrum to Far 
Ultraviolet Spectroscopic Explorer (FUSE) data taken towards 
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Table 3. Parity- and rotation-dependent inverse lifetimes for 
12 CO. 



Band 


A to t,f d 




Refs." 


# 


(s- 1 ) 


(s- 1 ) 




8 


3.6(1 1)+4.0(9)jc 


1.6(11)+13(10)jc 


1 


13 c 


3.4(10)+7.3(10)jc 




2 


16 


1.0(ll)+1.8(9)x 


1.0(1 l)+3.4(9)x 


1 


22 


1.83(9) 


1.91(9)+1.20(9)jc 


3 


25 


1.2(10) 


1.2(10)+2.4(9)x 


1 



a x sta nds for J'(J' + 1) . 

b (1) lEidelsberg et al I d2006l) : (2) lUbachs et alj il994l) ; (3) 

iDrabbels etal.U1993h . 
c This is a 'E + upper state, so there is no distinction between e and / 

parity. 



the star HD 203374A. Lastly, the Eid elsberg etajl (120041) value 
for the L'l ban d is 33 % l ower than that of | vDB88l but very sim- 
ilar to those of lER90l and lSheffer et al.1 d2003l) . so we adopt it as 
well. 

We compute the oscillator strengths of individual lines as 
the product of the appropriate Honl-London factor and the os- 
cillato r strength of the corresponding band dMorton & Noreaul 
Il994l) . Significant departures from Honl-London patterns have 
been reported fo r many N2 lines, someti mes even for the lowest 
rotational levels dStark et al.l2 005, 2008). The oscillator strength 
measured in one particular N2 band for the P(22) line was twenty 
times stronger than that for the P(2) line, due to strong mixing 
of the upper state with a nearby Rydberg state. For other bands 
where deviations from Honl-London factors were observed, the 
effect was generally less than 50% at J' = 10. Similar deviations 
are likely to occur for CO, but a lack of experimental data pre- 
vents us from including this in our model. Note, however, that 
large deviations are only expected for specific levels that happen 
to be strongly interacting with another state. Many hundreds of 
levels contribute to the photodissociation rate, so the effect of 
some erroneous individual line oscillator strengths is small. 

2.4. Lifetimes and predissociation probabilities 

Upon excitation, there is competition between dissociative and 
radiative decay. A band's predissociation probability (rf) can be 
computed if the upper st ate's to tal and radiative lifetimes are 
known: 77 = 1 — T>nt /t^h. IeR9Q| reported total lifetimes (l/A to t) 
for all their bands, but many of these are no more than order-of- 
magnitude estimates. Higher-resolution experime nts have since 



produced more accurate values for 17 of our bands ( Ubachs et al 



1994] 120051 ICacciani et al1fT998l I200TI l200l |Eidelsberg.etal 



2006). In several cases, values that differ by up to a factor of 
three are reported for different isotopologues. Where available, 
we take isotope-specific values. Otherwise, we follow the same 
procedure as for the oscillator strengths, and take the value of 
the isotopologue nearest in mass. 

The total lifetimes of some uppe r states have been 
shown to depend on the rotational level dDrabbels et al. 1993; 
lUbachs et al.1l2000t lEidelsberg et al.ll2006l) . In case of 'n states, 
a dependence on parity was sometimes observed as well. We in- 
clude these effects for the five bands in which they have been 
measured (Table[3]l. 

Recent experiments by Chakrabo rtv et al.l d2008l) suggest 
isotope-dependent photodissociation rates for the EQ, El, K0 
and W2 bands. These have been interpreted to imply different 



predissociation probabilities of individual lines of the various 
isotopologues due to a near-resonance accidental predissociation 
process. Similar effects have been reporte d for CIO2 and CO2 
ilLim etal.|[l999t iBhattacharva et al.ll2000T) . In this process, the 
bound-state levels into which the UV absorption takes place do 
not couple directly with the continuum of a dissociative state. 
Instead, they first transfer population to another bound state, 
whose levels happen to lie close in energy. For the CO E 1 state , 
this process was rotationally resolved by lUbachs et al.l ([2000) 
for all naturally occurring isotopologues and shown to be due 
to spin-orbit interaction with the k 3 Il v'—b state, which in turn 
couples with a repulsive state. The predissociation rates of the 
El state are found to increase significantly due to this process, 
but only for specific J' levels that accidentally overlap. For ex- 
ample, interaction occurs at J'-l, 9 and 12 for 12 CO, but at J'-\ 
and 6 for 13 C ls O. Since the dissociation probabilities for the El 
state due to direct predissociation were already high, 77 = 0.96, 
this increase in diss ociation rate only has a very minor effect 
(ICacciani et al.l 1998b . Moreover, under astrophysical conditions 
a range of J' values are pop ulated, so that th e effect of individ- 
ual levels is diluted. Since Chakrabortv et al.l did not derive line- 
by-line molecular parameters, we cannot easily incorporate their 
results into our model. In Sect. 16.31 we show that our results do 
not change significantly when we include the proposed effects in 
an ad-hoc way. 

The radiative lifetime of an excited state is a sum over the 
decay into all lower-lying levels, including the A l Yl and B '2 + 
electronic states and the v'^0 levels of the ground state. The 
total decay rate to the ground state is obtained by summing 
the o scillator strengths from Table Q] for each vibrational se- 
ries dMorton & Noreaulll994t ICacciani et al.|[l998l) . Theoretical 
work bv lKirbv & Cooper! (Il989l) showed that transitions to elec- 
tronic states other than the ground state contribute about 1% of 
the overall radiative decay rate for the C state and about 8% for 
the E state. No data are available for the excited states at higher 
energies. Fortunately, the radiative decay rate from these higher 
states to the ground state is small compared to the dissociation 
rate, so an uncertainty of ~10% will not affect the 77 values. 

The predissociati on prob abilities thus computed are nearly 
identical to those of vDB88: the largest difference is a 10% 
decrease for the CI band. This is due to the larger oscillator 
strength now adopted. 

There have been suggestions th at the CO state can also con- 
tribute to the photodissociation rate. ICacciani et al" I dlOOll) mea- 
sured upper-state lifetimes in the CO and CI states for several 
CO isotopologues. For the v'=0 state of I3 CO they found a total 
lifetime of 1770 ps, consistent with a value of 1780 ps for 12 CO, 
but different from the lifetime of 1500 ps in 13 C 18 0. Although 
the three values agree w ithin the mu tual uncertainties of 10-15% 
on each measurement. ICacciani et ail suggested that the heaviest 
species, 13 C ls O, has a predissociation yield of 77 = 0.17 rather 
than zero if the measurements are taken at face value and if the 
radiative lifetime of the CO state is presumed to be 1780 ps for all 
three sp ecies. The accurate a bsorption oscillator strength mea- 
sured bv lFederman et al.1 (1200 ll) for the CO band, 0.123 ± 0.016, 
implies a radiative lifetime that can be no longer than 1658 ps 
at the lower bound of measurement un certainty in / v /q . Taken 
together, the lifetime measu rements of |Cacciani et al.1 and the 
absorption measurements of iFederman et al.l favour aconserva- 
tive conclusion that the dissociation yield is zero for each of 
these three isotopologues. We assume the CO band is also non- 
dissociative in C 17 Q, C 18 Q and 13 C 17 Q; this is consistent with 
earlier studies (e.g., lvDB88llER90tlMo"rton & Noreaull 19941) . 
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2.5. Atomic and molecular hydrogen 

Lines of atomic and molecular hydrogen (H and H2) form an 
impor tant contribution to the overall shielding of CO. As in 
vDB88, we include H Lyman lines up to n-50 and H2 Lyman 
and Werner lines (transitions to the B l 2£ and C l Yl u states) 
from the v"=0, J" =0-1 levels of the electronic ground state. We 
adopt the line positio ns, oscillator strengths and lifetimes from 
Abgr all et alj dl993a|b|), as compiled for the freely available 
Meud on PDR code dLe Bourlot et al.|[l993l: iLe Petit et al.ll2002l 
2006)0 Ground-state rotational constants, required to compute 
the level populations, come from lJennings et ail d!984l) . 

3. Depth-dependent photodissociation 

3. 1 . Default model parameters 

The simplest way of modelling the depth-dependent photodisso- 
ciation involves dividing a one-dimensional model of an astro- 
physical object, irradiated only from one side, into small steps, 
in which the photodissociation rates can be assumed constant. 
We compute the abundances from the edge inwards, so that at 
each step we know the columns of CO, H, H2 and dust shielding 
the unattenuat ed radiation field. 

FollowinglvDB88jJLeBourlot et alJ (fl993l) . lLee et aD(fl996h 
and lLe Petit et alj (|2006), we treat the line and continuum atten- 
uation separately. For each of our 37 CO bands, we include all 
lines originating from the first ten rotational levels (/" =0-9) of 
the v"=0 level of the electronic ground state. That results in 855 
lines per isotopologue. In addition, we have 48 H lines and 444 
H2 lines, for a total of 5622. We use an adaptive wavelength 
grid that resolves all lines without wasting computational time 
on empty regions. For typical model parameters, the wavelength 
range from 91 1.75 to 1117.80 A is divided into -47,000 steps. 

We characterise the population distribution of CO over the 
rotational levels by a single temperature, r ex (CO). The H2 pop- 
ulation requires a more detailed treatme nt, because UV pump- 
ing pl ays a large role for the J">4 levels dvan Dishoeck & B lack 
1986). We populate the /"=0-3 levels according to a single tem- 
perature, r ex (H2), and adopt fixed columns of 4 x 10 15 , 1 x 10 15 , 
2 x 10 14 and 1 x 10 14 cirT 2 for J" =4-7. This reproduces ob- 
serv ed translucent cloud column d ensities to within a factor of 
two dvan Dishoeck & Blackl fl986 and references therein). The 
J">4 population scales with the UV intensity, so we re-evaluate 
this point for PDRs and circumstellar disks in Sects. I6~2l and l6. 31 

The line profiles of CO, H2 and H are taken to be Voigt func- 
tions, with default Dopp ler widths (b) of 0.3, 3.0 and 5.0 km s _1 , 
respectively. We adopt Drain^ d!978l) as our standard unattenu- 
ated interstellar radiation field. 

3.2. Unshielded photodissociation rates 

We obtain an unshielded CO photodisso ciation ra te of 2.6x 10~ 10 
s _1 . This rate is 30% higher than that of|vDB88, due to the gen- 
erally larger oscillator strengths in our data set. The new data for 
bands 33 and 24 (the E0 and K0 transitions) have the largest 
effect: they account for 63 and 21% of the overall increase. 
Clear ly, the rate dep e nds on the choice of radia tio n field. If we 
adoptlHabind (l968), GondhalekareTaD dl980t) or lMafhis et aT] 
( 1983) instead of Draine] d!978l) . the photodissociation rate be- 
comes 2.0, 2.0 or 2.3 x 10~ 10 s"\ respectively. lvDB88l reported 
the same relative differences between these three fields. 



2 http://aristote.obspm.fr/MIS/pdr/pdrl .html 



3.3. Shielding by CO, H 2 and H 

Self-shielding, shielding by H, H2 and the other CO isotopo- 
logues, and continuum shielding by dust all reduce the photodis- 
sociation rates inside a cloud or other environment relative to 
the unshielded rates. For a given combination of column depths 
(/V) and visual extinction (Ay), the photodissociation rate for iso- 
topologue i is 

h = Af&o,;®; exp(-yA v ) , (2) 

with x a scaling factor for the UV intensity and koj the unatten- 
uated rate in a given radiation field. The shielding function 0, 
accounts for self-shielding and shielding by H, H2 and the other 
CO isotopologues; tabulated values for some typical model pa- 
rameters are presented in Tables |5]-[8] The dust extinction term, 
exp(-yAy), is discussed in Sect. 13.41 Equation assumes the 
radiation is coming from all directions. If this is not the case, 
such as for a cloud irradiated only from one side, kyj should be 
reduced accordingly. 

For now, we ignore shielding by dust and compute the depth- 
dependent photodissociation rates due to line shielding only. 
Our test case is the centre of the diffuse cloud towards the star 
£ Oph. The observed column densities of H, H2, 12 CO and 
I3 CO are 5.2 x 10 20 , 4.2 x 10 20 , 2 5 x 10 15 and 1.5 x 10 13 
cirT 2 dvan Dishoeck & Blackl [l986t iLambert et alj 1 19941) . and 
we assume C n O, C ls O, 13 C n O and 13 C ls O column densities 
of 4.1 x 10 11 , 1.6 x 10 12 , 5.9 x 10 9 and 2.3 x 10 10 cirT 2 , consis- 
tent with observational constraints. For the centre of the cloud, 
we adopt half of thes e values. We set b (CO) = 0.48 km s _1 
and r ex (CO) = 4.2 K (ISheffer et al.ll 19921) . and we populate the 
H2 rotational levels explicitly according to the observed distri- 
bution. The cloud is illuminated by three times the Draine| field 
(X = 3). 

Table [4] lists the relative contribution of the most important 
bands at the edge and centre for each isotopologue, as well as the 
overall photodissociation rate at each point. The column densi- 
ties are small, but isotope-selective shielding already occurs: the 
I2 CO rate at the centre is lower than that of the other isotopo- 
logues by factors of three to six. 

The EQ band is the most im portant contributor at the edge. 
This was also found by vDB88, but the higher oscillator strength 
now adopted makes it even stronger. Going to the centre, it sat- 
urates rapidly for I2 CO: its absolute contribition to the total dis- 
sociation rate decreases by two orders of magnitude, and it goes 
from the strongest band to the 12th strongest b and. Th e five 
strongest bands at the centre are the same as in lvDB 88: Nos. 
13, 15, 19, 20 and 23. 

Figure [TJ illustrates the isotope-selective shielding. The left 
panel is centred on the R(l) line of the E0 band (No. 33 from 
Table [T). This line is fully saturated in 12 CO and the relative in- 
tensity of the radiation field (I/Io, with Iq the intensity at the 
edge of the cloud) goes to zero. 13 CO and C I8 also visibly re- 
duce the intensity, to I /Iq = 0.50 and 0.72, but the other three 
isotopologues are not abundant enough to do so. Consequently, 
these three are not self-shielded in the £ Oph cloud, but they 
are shielded by 12 CO, 13 CO and C 18 0. The weaker shielding of 
I3 C n O and 13 C ls O in the £0 band compared to C n O is due to 
their lines having less overlap with the 12 CO lines. 

The right panel of Fig. [TJ contains the R(0) line of the Wl 
band (No. 20), with the R(l) line present as a shoulder on the 
red wing. Also visible is the saturated B13 R(2) line of H2 at 
956.58 A. The Wl band is weaker than the E0 band, so 12 CO is 
the only isotopologue to cause any appreciable reduction in the 
radiation field and to be (partially) self-shielded. The shielding 
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Table 4. Relative and absolute shielding effects for the ten most important bands in the £ Oph cloud. 



Band 


33 


28 


24 


23 


22 


20 


19 


16 


15 


13 


Total c 


4) (A) b 


1076.1 


985.6 


970.4 


968.9 


968.3 


956.2 


950.0 


941.2 


940.0 


933.1 




12 CO 


Edge (%) d 


32.2 


4.7 


8.7 


3.0 


3.6 


3.4 


4.3 


5.1 


3.5 


3.1 


7.8(-10) 


Centre (%) d 


2.8 


0.4 


5.7 


9.3 


0.8 


12.1 


13.8 


6.1 


11.0 


7.0 


7.5(-ll) 


Shielding c 


0.0084 


0.0072 


0.063 


0.30 


0.022 


0.34 


0.31 


0.12 


0.30 


0.22 


0.10 


c 17 o 


Edge (%) 


32.0 


4.7 


8.7 


3.0 


3.6 


3.4 


4.3 


5.1 


3.5 


3.1 


7.8(-10) 


Centre (%) 


22.0 


0.1 


5.3 


3.5 


0.8 


8.7 


10.3 


14.2 


9.1 


2.4 


2.4(-10) 


Shielding 


0.21 


0.0067 


0.19 


0.36 


0.065 


0.79 


0.74 


0.85 


0.81 


0.23 


0.31 


c 18 o 


Edge (%) 


31.5 


5.1 


6.2 


1.5 


4.9 


3.7 


4.7 


5.2 


3.8 


3.4 


7.2(-10) 


Centre (%) 


35.3 


0.1 


4.7 


2.6 


7.8 


3.3 


7.7 


8.1 


5.9 


3.0 


3.7(-10) 


Shielding 


0.58 


0.0067 


0.39 


0.88 


0.82 


0.46 


0.83 


0.80 


0.81 


0.45 


0.52 


13 CO 


Edge (%) 


28.1 


6.5 


8.3 


1.8 


5.2 


3.4 


4.3 


4.9 


3.9 


2.4 


7.8(-10) 


Centre (%) 


18.5 


0.0 


5.3 


4.4 


11.2 


4.1 


10.7 


4.4 


9.1 


0.8 


2.8(-10) 


Shielding 


0.23 


0.0007 


0.23 


0.87 


0.76 


0.42 


0.87 


0.32 


0.83 


0.13 


0.35 


13 c 17 o 


Edge (%) 


28.0 


6.3 


8.3 


1.8 


5.2 


3.4 


4.3 


4.9 


3.9 


2.4 


7.8(-10) 


Centre (%) 


32.5 


0.0 


7.4 


2.9 


7.5 


5.0 


7.2 


8.5 


6.1 


1.0 


4.2(-10) 


Shielding 


0.63 


0.0007 


0.48 


0.88 


0.79 


0.78 


0.91 


0.94 


0.85 


0.24 


0.54 


13 C 18q 


Edge (%) 


30.1 


7.0 


5.9 


1.5 


4.7 


3.6 


4.4 


5.0 


4.0 


2.4 


7.5(-10) 


Centre (%) 


39.0 


0.0 


5.8 


2.4 


5.7 


1.7 


7.5 


8.5 


6.3 


1.7 


4.1(-10) 


Shielding 


0.71 


0.0007 


0.54 


0.87 


0.67 


0.26 


0.92 


0.94 


0.85 


0.38 


0.55 



a See the text for the adopted column densities, Doppler widths, excitation temperatures and radiation field. 
b 12 CO band head position. 

c Total photodissociation rate in s at the edge and the centre, and the shielding factor at the centre. 

d Relative contribution per band to the overall photodissociation rate at the edge and the centre of the cloud. 

e Shielding factor per band: the absolute contribution at the centre divided by the absolute contribution at the edge. 



relations. This causes the C I7 line to appear longwards of the 
C 18 line. 



3.4. Continuum shielding by dust 

Dust can provide a very strong attenuation of the radiation field. 
This effect is largely independent of wavelength for the 912— 
1 1 18 A radiation available to dissociate CO, so it affects all iso- 
topologues to the same extent. It can be expressed as an expo- 
nential function of the visual extinction, as expressed in Eq. (|2). 
For typical inter stellar dust grains (ra dius of 0. 1 fim and optical 
properties from lRoberge et all 1 19911), the e xtinction coefficient 
j is 3.53 for CO ( Ivan Dishoeck et alj 2006). Larger grains have 
less opacity in the UV and do not shield CO as strongly. For 
ice-coated grains w ith a mean radius of 1 fim, appropriate for 
circumstellar disks ( Jonkheid et al. 2006), the extinction coeffi- 
cient is only 0.6. The effects of dust shielding are discussed more 
fully in Sect. [6] 

Photodissociation of CO may still take place even in highly 
extincted regions. Cosmic rays or energetic electrons gener- 
ated by cosmic rays can excite H2, allowing it to emit in a 
multitude of bands, incl uding the Lyman and Werner systems 
dPrasad & Tara fdar 1983). The resulti ng UV photons can disso- 
ciate CO at a rate of about 10" 15 s" 1 dGredel et alJI 19871) . inde- 
pendent of depth. That is enough to increase the atomic C abun- 
dance by some three orders of magnitude compared to a situation 
where the photodissociation rate is absolutely zero. The cosmic- 
ray-induced photodissociation rate is sensitive to the spectro- 
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0.8 
0.6 
0.4 
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0.0 
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A-1075 (A) 



0.0 
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Fig. 1. Relative intensity of the radiation field (I/Iq) and intrinsic 
line profiles for the six CO isotopologues (<p, in arbitrary linear 
units) at the centre of the ( Oph cloud in two wavelength ranges. 



of the other five isotopologues is dominated by overlap with the 
12 CO and H2 lines. This figure also shows the need for accurate 
line positions: if the 13 C 18 line were shifted by 0.1 A in either 
direction, it would no longer overlap with the H2 line and be less 
strongly shielded. Note that the position of the Wl band has only 
been measured for 12 CO, 13 CO and C 18 0, so we have to compute 
the position for the other isotopologues from theoretical isotopic 
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scopic constants of CO, especially where it concerns the over- 
lap between CO an d H2 lines, s o it would be interesting to redo 
the calculations of iGredel et al.l with the new data from Table Q] 
However, that is beyond the scope of this paper. 

3.5. Uncertainties 

The uncertainties in the molecular data are echoed in the model 
results. When coupled to a chemical network, as in Sect. [6] the 
main observables produced by the model are the column den- 
sities of the CO isotopologues for a given astrophysical envi- 
ronment. The accuracy of the photodissociation rates is only 
relevant in a specific range of depths; in the average interstel- 
lar UV field, this range runs from an Ay of ~0.2 to ~2 mag. 
Photoprocesses are so dominant at lower extinctions and so slow 
at higher extinctions that the exact rate does not matter. In the in- 
termediate regime, both the absolute photodissociation rates and 
the differences between the rates for individual isotopologues are 
important. The oscillator strengths are the key variable in both 
cases and these are generally known rather accurately. Taking 
account of the experimental uncertainties in the band oscilla- 
tor strengths and of the theoretical uncertainties in computing 
the properties for individual lines, and identifying which bands 
are important contributors (Table [U, we estimate the absolute 
photodissociation rates to be accurate to about 20%. This error 
margin carries over into the absolute CO abundances and col- 
umn densities for the Ay w 0.2-2 mag range when the rates are 
put into a chemical model. The accuracy on the rates and abun- 
dances of the isotopologues relative to each other is estimated to 
be about 10% when summed over all states, even when we al- 
low for the k ind of isotope-spec i fic pre dissociation probabilities 
suggested bv lChakrabortv et all d2008l) . 

4. Excitation temperature and Doppler width 

The calculations of lvDB 88 were only done for low excita- 
tion temperatures of CO and H2. Here, we extend this work to 
higher temperatures, as required for PDRs and disks, and we re- 
examine the effect of the Doppler widths of CO, H2 and H on 
the photodissociation rates. We first treat four cases separately, 
increasing either r ex (CO), b(CO), T ex (H 2 ) or b(H 2 ). At the end 
of this section we combine these effects in a grid of excitation 
temperatures and Doppler widths. As a template model we take 
the centre of the £ Oph cloud, with column densities and other 
parameters as described in Sect. 13.31 

4.1. Increasing T ex (CO) 

As the excitation temperature of CO increases, additional rota- 
tional levels are populated and photodissociation is spread across 
more lines. Figure [2] visualises this for band 13 of 12 CO. At 4 K, 
only four lines are active: the R(0), R(l), P(l) and P(2) lines at 
933.02, 932.98, 933.09 and 933.12 A. The R(0) and P(l) lines 
are both fully self-shielded at the line centre. Going to 16 K, the 
R(0) line loses about 70% of its intrinsic intensity and ceases to 
be self-shielding. In addition, the R(2), P(3) and higher-/ lines 
start to absorb. The combination of less saturated low-/ lines and 
more active higher-/ lines yields a 39% higher photodissociation 
rate at 16 K compared to 4 K. 

A higher CO excitation temperature has the same favourable 
effect for I3 CO, which is partially self-shielded at the centre of 
the f Oph cloud. Its photodissociation rates increase by 16% 
when going from 4 to 16 K. C I8 is also partially self-shielded, 



but less so than 13 CO, so the favourable effect is smaller. At the 
same time, it suffers from increased overlap by 12 CO. The net 
result is a small increase in the photodissociation rate of 0.2%. 

The two heaviest isotopologues, 13 C n O and 13 C 18 0, are not 
abundant enough to be self-shielded. Their J"<2 lines generally 
have little overlap with the corresponding 12 CO lines, especially 
in the £0 band near 1076 A. This band, whose lines are amongst 
the narrowest in our data set, is the strongest contributor to the 
photodissociation rate at the centre of the cloud for 13 C n O and 
13 C ls O (Table @). In fact, its narrow lines are part of the rea- 
son it is the strongest contributor. The J" =3 and 4 lines that be- 
come active at 16 K do have some overlap with 12 CO. Without 
the favourable effect of less self-shielding, this causes the pho- 
todissociation rate for 13 C n O and 13 C 18 to decrease for higher 
excitation temperatures. The change is only small, though: 0.4% 
for 13 C n O and 2% for 13 C 18 0. 

Finally, C 17 experiences an increase of 18% in its pho- 
todissociation rate. Its lines lie closer to those of 12 CO than do 
the 13 C n O and 13 C ls O lines, so it is generally more strongly 
shielded. At 4 K, most of the shielding is due to the saturated 
R(0) lines of 12 CO. These become partially unsaturated at higher 
/ ex (CO), so the corresponding R(0) lines of C n O become a 
stronger contributor to the photodissociation rate, even though 
the shift towards higher-/ lines make them intrinsically weaker. 
Overall, increasing T ex (CO) from 4 to 16 K thus results in a 
higher C n O photodissociation rate. 



4.2. Increasing b(CO) 

The width of the absorption lines is due to Doppler broad- 
ening and natural (or lifetime) broadening. The integrated in- 
tensity in each line remains the same when b(CO) increases, 
so a larger width is accompanied by a lower peak intensity. 
The resulting reduction in self-shielding then causes a higher 
I2 CO photodissociation rate, as shown in Fig. [2] for band 13. 
However, the effect is rather small because the Doppler width 
is smaller than the natural width for most lines at typical b val- 
ues. Natural broadening is the dominant broadening mechanism 
up to b(CO) « 6 x 10~ 12 A tot , with both parameters in their nor- 
mal units. The R(0) line of band 13 has an inverse lifetime of 
1.8 x 10 11 s _1 (TablesQ]and[3]), so Doppler broadening becomes 
important at about 1 km s . From 0.3 to 3 km s _1 , as in Fig. [2] 
the line width only increases by a factor of 1.9. Integrated over 
all lines, the 12 CO photodissociation rate becomes 26% higher. 

The rates of the other five isotopologues decrease along this 
b(CO) interval due to increased shielding by the £0 lines of 
I2 CO. With an inverse lifetime of only 1 x 10 9 s _1 , Doppler 
broadening is this band's dominant broadening mechanism in 
the regime of interest. A tenfold increase in the Doppler param- 
eter from 0.3 to 3 km s _1 results in a nearly tenfold increase 
in the line widths. At 0.3 km s _1 , the £0 lines of 12 CO are 
still sufficiently narrow that they do not strongly shield the lines 
of the other isotopologues. This is no longer the case at 3 km 
s _1 . I3 CO still benefits somewhat from reduced self-shielding 
in other bands, but it is not enough to overcome the reduced 
strength of the £0 band, and its photodissociation rates decrease 
by 2%. The decrease is 13% for C I7 and 26-28% for the re- 
maining three isotopologues. The relatively small decrease for 
C n O is due to its £0 band being already partially shielded by 
I2 CO at 0.3 km s _1 , so the stronger shielding at 3 km s _1 has less 
of an effect. 
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Fig. 2. Illustration of the effect of increasing r ex (CO), b(CO), r ex (H2) or b(}\2) (from left to right) in our ( Oph cloud model (Sect. 
13 - 31 > - Top: relative intensity of the radiation field, including absorption by only 12 CO (left two panels) or H2 (right two panels). 
Middle: intrinsic line profile for band 13 (933.1 A) or 30 (1002.6 A) of 12 CO. Bottom: photodissociation rate per unit wavelength, 
multiplied by a constant as indicated. 



4.3. Increasing T ex (U 2 ) orb(H 2 ) 



Increasing the excitation temperature of H2, while keeping the 
CO parameters constant, results in a decreased photodissocia- 
tion rate for all six isotopologues. The cause, as again illustrated 
in Fig. |2 is the activation of more H2 lines. At r ex (H2) = 10 
K, the R(l) line of the B8 band at 1002.45 A is very narrow 
and does not shield the F0 band (No. 30) of the CO isotopo- 
logues. (The continuum-like shielding visible in Fig. [2]is due to 
the strongly saturated B8 R(0) line at 1001.82 A.) It becomes 
much more intense at 30 K and widens due to being saturated, 
thereby shielding part of the FQ band. The same thing happens 
to other CO bands, resulting in an overall rate decrease of 0.6- 
2.5%. There is no particular trend visible amongst the isotopo- 
logues; the magnitude of the rate change depends purely on the 
chance that a given CO band overlaps with an H2 line. 

Similar decreases of one or two percent in the CO photodis- 
sociation rates are seen when the H2 Doppler width is changed 
from 1 to 10 km s _1 . As the H2 lines become broader, the amount 
of overlap with CO increases across the entire wavelength range. 
As an example, Fig.|2]shows again the 1002.0-1002.8 A region, 
where the B8 R(l) line of H2 further reduces the contribution of 
the F0 band to the 12 CO photodissociation rate. 



4.4. Grid of T ex and b 

We now combine the four individual cases into a grid of exci- 
tation temperatures and Doppler widths to see how they influ- 
ence each other. r ex (CO) is raised from 4 to 512 K in steps of 
factors of two. The v"=l vibrational level of 12 CO lies at 2143 
crrr 1 above the v"-Q level, so it starts to be thermally popu- 
lated at ~500 K. No data are available on dissociative transi- 
tions out of this level, so we choose not to go to higher exci- 
tation temperatures. We increase the range of rotational levels 
up to /"=39, at 2984 cirT 1 above the J"=0 level for I2 CO. At 
r ex (CO) = 512 K, the normalised population distribution peaks 
at J"=9 and decreases to 9.7 x 10~ 5 at /"=39. The H 2 excita- 
tion temperature is set to [T ex (CO)] to take account of the fact 
that its critical densities for thermalisation are lower than those 
of CO. Where necessary, absorption by rotational levels above 
our normal limit of J" -1 and by non-zero vibrational le vels is 
taken into account dDabrowskill984tlAbgrall et al.ll993allbh . All 
H2 rovibrational levels are strictly thermally populated; no UV 
pumping is included. The grid is run for CO Doppler widths of 
0.1, 0.3, 1.0 and 3.2 km s -1 ; we set b(H 2 ) = Vl4Z>(CO) and 
b(H) = V28b(CO), corresponding to the differences appropriate 
for thermal broadening. 

The left set of panels in Fig. [3] shows the photodissociation 
rate of the six isotopologues at the centre of the £ Oph cloud 
as a function of excitation temperature for the different Doppler 
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Fig. 3. Left: photodissociation rate of the CO isotopologues as a function of excitation temperature, normalised to the rate at 4 K, 
at four different Doppler widths, for parameters corresponding to the centre of the f Oph diffuse cloud. Right: photodissociation 
rates as a function of Doppler width, normalised to the rate at 0. 1 km s _I , at four different excitation temperatures. The other model 
parameters are as in Sect. [33] except b(H 2 ) = Vl4Z?(CO), b(H) = V28£(CO) and r ex (H 2 ) = [r ex (CO)] L5 . 




widths. The rates are normalised to the rate at 4 K. The CO 
rate increases from 4 to 16 K, as described in Sect. 14.11 At 
higher temperatures the increased overlap with H2 lines takes 
over and the rate goes down. As long as the CO excitation tem- 
perature is less than ~100 K, the 12 CO rate remains constant up 
to b(CO) = 0.3 km s _1 and increases as b(CO) goes from 0.3 to 
3.2 km s _I (right side of Fig. [3]). At higher temperatures there is 
so much shielding by H2 that reduced self-shielding in the CO 
lines has no discernable effect on the rate. Instead, the rate goes 
down with b due to stronger shielding by the broadened H2 lines. 



The CO rate also increases initially with r ex (CO) and then 
goes down as H2 shielding takes over. The rate increases from 
b(CO) — 0.1 to 0.3 km s , but decreases for higher values as 
described in Sect. 14.21 For the remaining four isotopologues, 
the plotted curves likewise result from a combination of weaker 
shielding by 12 CO and stronger shielding by H2. At CO excita- 
tion temperatures between 4 and 8 K, the rates typically change 
by a few per cent either way. Going to higher temperatures, all 
rates decrease monotonically. Likewise, the rates generally de- 
crease towards higher b(CO) values. 
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A change in behaviour is seen when increasing r ex (CO) from 
256 to 5 1 2 K. It is at this point that the v" >0 levels of H2 become 
populated. Less energy is now needed to excite H2 to the B and 
C states, so absorption shifts towards longer wavelengths. This 
causes even stronger shielding in the heavy CO isotopologues, 
for whom the E0 band at 1076 A is still an important contributor 
to the photodissociation rate, at least as long as b(CO) does not 
exceed 0.3 km s" 1 . The E0 band is strongly self-shielded in 12 CO 
(Table|4]i, so the shift of the H2 absorption to longer wavelengths 
does not reduce its contribution by much. In fact, the weaker H2 
absorption at shorter wavelengths allows for an increased contri- 
bution of bands like Nos. 13 and 16 at 933 and 941 A, causing 
a net increase in the 12 CO photodissociation rate from 256 to 
512 K. The situation changes somewhat when the CO Doppler 
width increases to 1 .0 km s" 1 or more. The E0 band of the heavy 
isotopologues is now much less of a contributor, because it is 
shielded by the broader 12 CO lines. The shifting H2 absorption 
does not cause any additional shielding, so the rates remain al- 
most the same. Furthermore, the H2 lines are also broader and 
continue to shield the I2 CO bands at shorter wavelengths, pre- 
venting its photodissociation rate from increasing like it does in 
the low-b cases. 

The main astrophysical consequence becomes clear when 
we look at the photodissociation rates of the five heavier iso- 
topologues with respect to that of 12 CO. Figure |4] shows this 
ki/k( l2 CO) ratio as a function of r ex (CO) and b(CO). In these 
plots, the I2 CO rate shows as a horizontal line with a value 
of unity. The ratios generally decrease with both parameters: a 
higher excitation temperature and a larger Doppler width both 
cause less self-shielding in 12 CO, so the rate differences between 
the isotopologues become smaller. Shielding by H2 increases at 
the same time, further reducing the differences between the iso- 
topologues. This means that photodissociation of CO is more 
strongly isotope-selective in cold sources than in hot sources. 

5. Shielding function approximations 

It is unpractical for many astrophysical applications to do the full 
integration of all 5622 lines in our model every time a photodis- 
sociation rate is required. Therefore, we present approximations 
to the shielding functions introduced in Eq. ©. The approxi- 
mations are derived for several sets of model parameters and are 
valid across a wide range of astrophysical environments (Sect. 

E2). 

5.1. Shielding functions on a two-dimensional grid ofN(CO) 
andN(H 2 ) 

The transition from atomic to molecular hydrogen occurs much 
closer to the edge of the cloud than the C + -C-CO transition, so 
the column density of atomic H is roughly constant at the depths 
where shielding of CO is important. In addition, H shields CO 
by only a few per cent. Therefore, it is a good approximation to 
compute the shielding functions on a grid of CO and H2 column 
densities, while taking a constant column of H. It is sufficient 
to express the shielding of all CO isotopologues as a function 
of 7V( 12 CO), because self-shielding of the heavier CO isotopo- 
logues is a small effect compared to shielding by 12 CO. 

T able [5] presents the shielding functions in the same man- 
ner as vDB88 did, but for somewhat different model parameters: 
b(CO) = 0.3 instead of 1.0 km s _1 , 7/ ex (CO) = 5 instead of 
10 K, and r ex (H 2 ) = 5 1 5 instead of 10 15 K. The present pa- 
rameters correspond more closely to what is observed in dif- 



fuse and translucent clouds. The column density ratios for the 
six isoto pologues are ke pt constant at the elemental isotope ra- 
tios from IWUsoiil ^199§: [ 12 C]/[ 13 C] = 69, [ 16 0]/[ 18 0] = 557 
and [ 18 0]/[ 17 0] = 3.6. A small column of 5.2 X 10 15 cm 2 of 
H 2 at J" =4-7 is included throughout (except at log/V(H 2 ) = 0) 
to account for UV pumping (Sect. [XTT ipl Shielding functions for 
larger values of b and T ex and for other isotope ratios are given 
in Tables |6]-[8] in the online appendix to this paper. For ease of 
use, we have also set up a webpag^] where the shielding func- 
tions can be downloaded in plain text format. This webpage of- 
fers shielding functions for a wider variety of parameters than 
is possible to include in this paper. In addition, it uses a grid of 
7V( 12 CO) and N(H2) values that is five times finer than the grid 
in Tables |5]-[8] allowing for more accurate interpolation. 

For column densities of up to 10 17 cm" 2 of CO and 10 21 cm -2 
of H2, our shie lding functions are generally within a few per cent 
of the vDB88 values when corrected for the difference in b and 
r ex . Larger differences occur for larger columns: we predict the 
shielding to be about five times weaker at N{CO) = 10 19 cm -2 
and more than a hundred times stronger at N(H2) = 10 23 cm" 2 . 
The 912 -1118 A wavelength range was divided into 23 bins by 
vDB88, and most lines were included only in one bin to speed 
up the computation. We integrate all lines over all wavelengths. 
As the lines get strongly saturated at high column depths, ab- 
sorption in the line wings becomes important. Thus, H2 lines can 
cause substantial shielding over a range of more than 10 A, while 
CO lines may still absorb several A away from the line centre. 
The binned integration method of vDB88 did not take these ef- 
fects into account, so they underpredicted shielding at large H2 
columns and overpredicted shielding at large CO columns (Fig. 
[5). It should be noted, however, that photodissociation at these 
depths is typically already so slow a process that it is no longer 
the dominant destruction pathway for CO. In addition, a large 
CO column is usually accompanied by a large H2 column, so the 
two effects partially cancel each other. 

5.2. Comparison between the full model and the 
approximations 

Despite being computed for a limited number of model param- 
eters, the shielding functions from Tables [5}{8] provide a good 
approximation to the rates from the full model for a wide range 
of astrophysical environments. Section 16.11 presents a grid of 
translucent cloud models, where the photodissociation model is 
coupled to a chemical network and CO is traced as a function of 
depth. This presents a large range of column densities, with the 
ratios between the isotopologues deviating from the fixed val- 
ues adopted for Tables [5|-[8] The grid covers gas densities from 
100 to 1000 cm" 3 and gas temperatures from 15 to 100 K, while 
keeping the excitation temperatures and Doppler widths constant 
at the values used for Table|5] Altogether, the grid contains 2880 
points per isotopologue where the photodissociation rate is com- 
puted. Here, we compare the photodissociaton rates from the ap- 
proximate method to the full integration for each of these points. 
The sensitivity of k( l3 CO) to the 7V( 12 CO)//V( 13 CO) ratio can be 
corrected for in a simple manner: we use the shielding functions 
from Table|8]when the ratio is closer to 35 than to 69, and those 
from Table |5]otherwise. 

The rate from our approximate method is within 10% of the 
"real" rate in 98.3% of all points (Fig. [6}. In no cases is the dif- 

3 Although not mentioned explicitly by vDB8j§, their tabulated 
shielding functions also include this extra column of J" =4-7 H 2 . 

4 http://www.strw.leidenuniv.nl/~ewine/photo 



R. Visser et al.: The photodissociation and chemistry of interstellar CO isotopologues 



11 



Table 5. Two-dimensional shielding functions @[N( l2 CO),N(U 2 )] for b(CO) = 0.3 km s -1 , r ex (CO) = 5 K and 
7V( 12 CO)//V( 13 CO) = 69. a 
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12 CO: unattenuated rate koj = 


2.592 x 10- 


10 g-l 









1.000 


8.080 


-i) 


5.250(-l) 


2.434(-l) 


5.467(-2) 


1.362(-2) 


3.378(-3) 


5.240(-4) 


19 


8.176(-1) 


6.347 


-1) 


3.891(-1) 


1.787(-1) 


4.297(-2) 


1.152(-2) 


2.922(-3) 


4.662(-4) 


20 


7.223(-l) 


5.624 


-1) 


3.434(-l) 


1.540(-1) 


3.515(-2) 


9.23l(-3) 


2.388(-3) 


^i o r\r\ / a \ 

3.899(-4) 


21 


3.260(-l) 


2.810 


-1) 


1.953(-1) 


8.726(-2) 


1.907(-2) 


4.768(-3) 


1.150(-3) 


1.941(-4) 


22 


1.108(-2) 


1.081 


-2) 


9.033(-3) 


4.441 (-3) 


1.102(-3) 


O C A A t A\ 

2.644(-4) 


7.329(-5) 


1.437(-5) 


23 


3 938f-71 


3.938 


-7) 


3.936(-7) 


3.923(-7) 


3.901(-7) 


3 893( , -7 , l 


3.890(-7) 


3 875f-7~) 










3: unattenuated rate koj = 


2.607 x 10- 


I0 -l 

s 









1.000 


9.823 


-i) 


8.911(-1) 


6.149(-1) 


3.924(-l) 


2.169(-1) 


4.167(-2) 


2.150(-3) 


19 


8.459(-l) 


8.298 


-1) 


7.490(-l) 


5.009(-l) 


3.196(-1) 


l.850(-l) 


3.509(-2) 


1.984(-3) 


20 


7.337(-l) 


7.195 


-1) 


6.481(-1) 


4.306(-l) 


2.741(-1) 


l.556(-l) 


2.645(-2) 


1.41 1(-3) 


21 


3.335(-l) 


3.290 


-1) 


3.039(-l) 


2.293(-l) 


1.685(-1) 


9.464(-2) 


1.460(-2) 


6.823(-4) 


22 


1.193(-2) 


i i n i 


-2) 


1.172(-2) 


1.095(-2) 


9.395(-3) 


C £.A A I 1\ 

J.644(-3) 


1.183(-3J 


2.8JD(-D) 


23 


3.959(-7) 


3.959 


-7) 


3.959(-7) 


3.959(-7) 


3.958(-7) 


3.954(-7) 


3.924(-7) 


3.873(-7) 








C 18 


3: unattenuated rate koj = 


2.392 x 10- 


io s -i 









1.000 


9.974 


-1) 


9.777(-l) 


8.519(-1) 


5.060(-l) 


l.959(-l) 


2.764(-2) 


1.742(-3) 


19 


8.571(-1) 


8.547 


-1) 


8.368(-l) 


7.219(-1) 


4.095(-l) 


l.58l(-l) 


2.224(-2) 


1.618(-3) 


20 


7.554(-l) 


7.532 


-1) 


7.371(-1) 


6.336(-l) 


3.572(-l) 


1.372(-1) 


1.889(-2) 


1.383(-3) 


21 


3.559(-l) 


3.549 


-1) 


3.477(-l) 


3.035(-l) 


1.948(-1) 


7.701(-2) 


1.071(-2) 


6.863(-4) 


22 


1.214(-2) 


1.212 


-2) 


1.199(-2) 


1.105(-2) 


8.233(-3) 


3.324(-3) 


6.148(-4) 


3.225(-5) 


23 


4.251(-7) 


4.251 


-7) 


4.251(-7) 


4.251 (-7) 


4.249(-7) 


4.233(-7) 


4.180(-7) 


4.142(-7) 



CO: unattenuated rate koj 






1.000 


9.824( 


-1) 


9.019(-1) 


6.462(-l) 


3.547(-l) 


9.907(-2) 


1.131(-2) 


7.591(-4) 


19 


8.447(-l) 


8.276( 


-1) 


7.502(-l) 


5.113(-1) 


2.745(-l) 


7.652(-2) 


8.635(-3) 


6.747(-4) 


20 


7.415(-1) 


7.266( 


-1) 


6.581(-1) 


4.45 1(-1) 


2.360(-l) 


6.574(-2) 


7.187(-3) 


5.429(-4) 


21 


3.546(-l) 


3.502( 


-1) 


3.270(-l) 


2.452(-l) 


1.398(-1) 


3.750(-2) 


3.973(-3) 


2.703(-4) 


22 


1.180(-2) 


1.177( 


-2) 


1.153(-2) 


1.023(-2) 


6.728(-3) 


1.955(-3) 


2.665(-4) 


1.471(-5) 


23 


2.385(-7) 


2.385( 


-7) 


2.385(-7) 


2.384(-7) 


2.379(-7) 


2.348(-7) 


2.310(-7) 


2.292(-7) 








J C 17 


O: unattenuated rate k i 


= 2.598 x 10- 


10 s -i 









1.000 


9.979( 


-1) 


9.820(-l) 


8.832(-l) 


5.942(-l) 


3.177(-1) 


1.523(-1) 


3.885(-2) 


19 


8.540(-l) 


8.520( 


-1) 


8.374(-l) 


7.469(-l) 


4.901(-1) 


2.677(-l) 


1.302(-1) 


3.135(-2) 


20 


7.405(-l) 


7.387( 


-1) 


7.254(-l) 


6.439(-l) 


4.198(-1) 


2.333(-l) 


1.142(-1) 


2.607(-2) 


21 


3.502(-l) 


3.494( 


-1) 


3.434(-l) 


3.076(-l) 


2.214(-1) 


1.386(-1) 


6.941(-2) 


1.195(-2) 


22 


1.279(-2) 


1.278( 


-2) 


1.267(-2) 


1.198(-2) 


1.045(-2) 


7.743(-3) 


4.088(-3) 


4.581(-4) 


23 


2.370(-7) 


2.370( 


-7) 


2.370(-7) 


2.370(-7) 


2.369(-7) 


2.368(-7) 


2.359(-7) 


2.312(-7) 








3 C 18 


O: unattenuated rate koj 


= 2.503 x 10" 


io s -i 









1.000 


9.988( 


-1) 


9.900(-l) 


9.329(-l) 


7.253(-l) 


3.856(-l) 


1.524(-1) 


2.664(-2) 


19 


8.744(-l) 


8.734( 


-1) 


8.656(-l) 


8.164(-1) 


6.403(-l) 


3.441(-1) 


1.347(-1) 


2.491(-2) 


20 


7.572(-l) 


7.562( 


-1) 


7.492(-l) 


7.047(-l) 


5.518(-1) 


3.006(-l) 


1.185(-1) 


2.224(-2) 


21 


3.546(-l) 


3.542( 


-1) 


3.506(-l) 


3.283(-l) 


2.638(-l) 


1.666(-1) 


6.887(-2) 


1.149(-2) 


22 


1.561(-2) 


1.560( 


-2) 


1.550(-2) 


1.475(-2) 


1.235(-2) 


7.850(-3) 


3.416(-3) 


5.290(-4) 


23 


2.490(-7) 


2.490( 


-7) 


2.490(-7) 


2.489(-7) 


2.487(-7) 


2.482(-7) 


2.471 (-7) 


2.421(-7) 



2.595 x 10- 



These shielding functions were computed for the lDraing d!978h radiation field (x = 1) and the following set of parameters: b(CO) = 0.3 km 
s-\ fc(H 2 ) = 3.0 km s- 1 and fc(H) = 5.0 km s _1 ; r ex (CO) = 5 K and T CX (H 2 ) = 5 1 5 K; 2V(H) = 5 x 10 20 cm 4 ; JV( l2 C0)/iV( l3 C0) = 69, 
A^( 12 CO)/Af(C 18 0) = Af( 13 CO)/iV( 13 C 18 0) = 557 and Af(C 18 0)MC n O) = M 13 C ls O)/Af( 13 C 17 0) = 3.6. Self-shielding is mostly negligible 
for the heavier isotopologues, so all shielding functions are expressed as a function of the 12 CO column density. Continuum attenuation by 
dust is not included in this table (see Eq. (O). 



ference between the approximate rates and the full model more 
than 40%. Perhaps even more important than the absolute pho- 
todissociation rates are the ratios, Q, between the rates of 12 CO 
and the other five isotopologues ^CO: 

Q ~ (3) 

The shielding functions from Tables|5]and[8]together reproduce 
the ratios from the full method to the same accuracy as the abso- 
lute rates. At 98.4% of all points in the grid of translucent clouds, 
the ratios are off by less than 10% (Fig. [6]). Similar scores can be 



obtained for models of PDRs or other environments if one uses 
shielding functions computed for the right combination of pa- 
rameters. Otherwise, the accuracy goes down. For example, the 
shielding functions from Table [5] (r ex (CO) = 5 K) can easily 
give photodissociation rates off by a factor of two when applied 
to a high-density, high-temperature PDR. 

6. Chemistry of CO: astrophysical implications 

Photodissociation is an important destruction mechanism for CO 
in many environments. In this section, we couple the photodis- 
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Fig. 5. Illustration of the effect of not integrating strongly satu- 
rated lines over the entire wavelength regime. Full black curve: 
relative intensity of the radiation field at N(CO) = 10 18 and 
NQI2) = 10 22 cirT 2 for the parameters of Table[5] computed with 
a full integration of all lines over all wavelengths. Red curve: 
the sam e, but computed with the binned integration method of 
lvDB88l Dotted curve: attenuation due to the C3 R (0) line o f H 2 , 
centred on 946.42 A. This line was not included by vDB88 when 
computing the shielding in the 938.9-943.1 A bin, thus overes- 
timating the photodissociation rate in the W2 band at 941 . 1 A. 



abs(Q 
10 




Fig. 6. Histogram of the absolute relative difference between the 
approximate method (Tables [5] and [8j see text for details) and 
the full computation for the absolute photodissociation rates (k, 
black line) and the isotopologue rate ratios (Q, red line) in a grid 
of translucent cloud models. The data for all six isotopologues 
are taken together, but the Q values for 12 CO (which are unity 
regardless of the method) are omitted. 



sociation model to a small chemical network in order to ex- 
plore abundances and column densities. Specifically, we model 
the CO chemistry in translucent clouds, PDRs and circumstel- 
lar disks, and we compare our results to observations of such 
objects. The photodissociation rates are computed with our full 
model throughout this section. 

6.1. Translucent clouds 
6.1.1. Model setup 

Translucent clouds, with visual extinctions between 1 and 5 mag, 
form an excellent test case for our CO photodissociation model. 



Two recent studies of several dozen lines of sight through diffuse 
(Ay < 1 mag) and translucent clou ds provide a set of CO and 
H2 column densitie s for comparison dSonnentrucker et ail 2007; 
Sheff er et alJl2008T) . These surveys show a clear correlation be- 
tween ./V(H2) and A^( I2 CO), with distinctly differ ent slopes for 
H2 columns of less and more than 2.5 x 10 20 cirr 2 . [Sheffer et all 
attributed this break to a change in the formation mechanism of 
CO. Their models show that the two-step conversion from C + to 
CH + and CO + , 



H 7 



CH + 



H 



CH + 



O 



CO + + H 



(4) 
(5) 



followed by reaction with H (forming CO directly) or H2 (form- 
ing HCO + , which then recombines with an electron to give CO) 
is the dominant pathway at low column densities. However, the 
highly endothermic Reaction (0]i is not fast enough at gas ki- 
netic temperatures typical for these environments to explain the 
observed abundances of CH + and CO. Suprather mal chemistry 
has been suggested as a s olution to this problem . ISheffer et alj 
followed the approach of Feder man et al.1 (119961) . who argued 
that Alfven waves entering the cloud from the outside re- 
sult in non-thermal motions between ions and neu t rals. O ther 
mechanisms have been suggested by iJoulain et al. (1998]) and 
iPetv & Falgarone (2000). The effect of the Alfven waves can 
be incorporated into a chemical model by replacing the kinetic 
temperature in the rate equation for Reaction (|4]i and all other 
ion-neutral reactions by an effective temperature: 



3k K 



(6) 



Here, k^ is the Boltzmann constant, fi is the reduced mass of the 
reactants and va is the Alfven speed. The Alfven waves reach a 
depth of a few 10 20 crrT 2 of H2, corresponding to an A v of a few 
tenths of a magnitude, beyond which suprathermal chemistry 
ceases to be important. CO can therefore no longer be formed 
efficiently through Reactions © and (|5), and the reaction be- 
tween C + and OH (producing CO either directly or via a CO + 
intermediate) takes over as the key route to CO. The identifica- 
tion of the se two different chemistry r egimes supports the con- 
clusion of Zsargo & Federmanl (120031) that suprathermal chem- 
istry is required to explain observed CO abundances in diffuse 
environments. Suprathermal chemistry al so drives up the HCO + 
abun dance, con firming the conclusion of iLiszt & Luc as (1994) 
and[Liszt (2007) that HCO + is the dominant precursor to CO in 
diffuse clouds. 

We present here a grid of translucent cloud models to see 
how well the new photodissociation results match the observa- 
tions. We set the Alfven speed to 3.3 km s _1 for N( U 2 ) < 4xl0 2Q 
cm" 2 and to zero for larger column densities dSheffer et al.l 
2008). The grid comprises densities («h = «(H) + 2«(H2)) of 
100, 300, 500, 700 and 1000 cm 4 , gas temperatures of 15, 30, 
50 and 100 K, and relative UV intensities (x) of 1, 3 and 10. The 
dust temperature is assumed to stay low for all models (~ 15 K), 
so the H2 formation rate does not change. The ionisation rate of 
H due to cosmic rays is set to a constant value of 1.3 x 10" 17 
s . Attenuation by 0.1 /im dust grains (Sect. [3~4l is taken into 
account. The Doppler widths and level populations are as de- 
scribed in Sect. 13.11 with CO and H2 excitation temperatures of 
5 and 5 15 K. Taking other b or T ex values plausible for these 
environments does not alter our results significantly. All models 
are run to an Ay of 5 mag; results are also presented for a range 
of smaller extinctions. 



R. Visser et al.: The photodissociation and chemistry of interstellar CO isotopologues 



13 



Table 9. Elemental abundances. 



Element 


Abundance relative to «h 


He 


1.00 x 10-' 


12 C 


1.40 x 10 4 


13 C 


2.03 x 10- 6 




3.10 x 10" 4 


17 o 


1.55 x 10- 7 


18Q 


5.57 x 10- 7 



The models require a chemical network to compute the abun- 
dances at each depth step. Since we are only interested in CO, 
the number of relevant species and reactions is limit ed. We adopt 
the ne twork from a recent PDR benchmark study dRollig et alj 
|2007|) . which includes only 3 1 species consisting of H, He, C 
and O. We duplicate all C- and O-containing species and re- 
actions for 13 C, n O and ls O. Freeze-out and thermal evapo- 
ration are added for all neutral species, but no grain-surface 
reactions are included other than H2 formation according to 
Blac k & van Dishoeckl dl987l) . We add ion-molecule exchange 
reactions such as 



12 co + 13 c + 



Z C + + u CO + 35 K, 



H 12 CO + + 13 CO I2 CO + H 13 CO + + 9K, 



(7) 
(8) 



which can enhance the abundances of the heavy isotopologues 
of CO and H CO + (|mtsgjietaD[i97l ISmith & Adamsl[T980t 
lLanger etall Il984l) . The temperatur e-dep e ndenc e of the rate 
of these two reacti ons was fitted by|L iszt (2007); the alterna- 
tive equations from IWoods & Willacvl d2009l) give the same re- 
sults. The effective temperature from Eq. © is used instead 
of the kinetic temperature for all ion-neutral reactions, includ- 
ing Reactions (0 and dB). Altogether, the network contains 
118 speci es and 1723 re a ctions . We adopt the elemental abun- 
dances of ICardelli et alj d 19961) and isotope ratios appropriate 
for the local ISM ([ 12 C1/[ 13 C1 = 69, [ 16 0]/[ 18 0] = 557 and 
[ 18 0]/[ 17 0] = 3.6;;Wilson| fT999l) : the complete list of elemental 



abundances is given in Table[9] Chemical steady state is reached 
at all depths after ~ 1 Myr, regardless of whether the gas starts in 
atomic or molecular form. 

6.1.2. 12 CO 

The left panel of Fig. shows N( 12 CO) versus N(H 2 ) for all 
depth steps in our grid. These and other column densities are also 
listed for a set of selected Ay values in Tables [TOlFOl in the online 
appendix to this paper. The scatter in the data is due to the dif- 
ferent physical parameters. For a given N(H2), iV( 12 CO) is about 
two times larger at r gas = 100 K than at 15 K. The formation rate 
of H2 increases with temperature; through the chain of reactions 
starting with Reaction ©, that results in a larger CO abundance 
and column density. Up to NQI2) = 10 21 cm -2 , increasing the 
gas density by a factor of ten increases N( l2 CO)/N(H2) also by 
about a factor of ten. This is due to the photodissociation rate 
being mostly independent of density, while the rates of the two- 
body reactions forming CO are not. Photodissociation ceases to 
be the main destruction mechanism for CO deeper into the cloud, 
so increasing «e has a smaller effect there. Increasing the UV 
intensity from^ = 1 to 10 has the simple effect of decreasing 
N( n CO)/N(U 2 ) roughly tenfold for N(U 2 ) < 10 21 cnr 2 . For 
larger depths, changing x only has a small effect. These depen- 
dencies o n the physical param eters are consistent with the obser- 
vations of Sheff er et al J dl008). 



The full set of points already gives some indication that the 
data are not well represented by a global single-slope correla- 
tion. The right panel of Fig. [7] shows the 100-point means of 
th e full set, rev e aling the same two distinct regimes as found 
bv lShefferetai] (2008). Taking the uncertainties in both the ob- 
servations and the models into account, the two sets of power- 
law fits (indicated by the black and grey lines) are identical. The 
break between the two slopes occurs at N{\\j) = 2.5 x 10 20 cm -2 
and is due to the switch from suprathermal to normal chem- 
istry. It is clearly unrelated to the switch from UV photons to 
He + or as the main destroyer of 12 CO, which does not oc- 
cur until N(\\i) = 1-2 x 10 21 cm 2 . When we extend the grid 
to higher Ay (not plotted), we quickly reach the point where the 
12 CO abundance equals the elemental 12 C abundance and the 
plot of Af( 12 CO) versus MH2) continues with a slope of unity. 
The model results from vDB88 (dashed lines, corresponding to 
X = 0.5, 1 and 10) do not extend to low enough N(H2) to show 
the suprathermal regime, but they do show at the upper end the 
transition to the regime where all gas-phase carbon is in CO. 



6.1.3. 13 CO 

With our model able to reproduce the observational relationship 
between H2 and 12 CO, the next step is to look at the heavie r iso- 
topologues.lSonnentrucker et alj d2007l) . lBurgh et all d2007l) and 
Shef fer et alJ d2007l) together presented a sample of 29 diffuse 
and translucent sources with derived column densities for both 
12 CO and 13 CO. These data are plotted in the left panel of Fig. 
EKgrey crosses) as N( n CO) versus A^Hb). The 100-point means 
from our model grid (black squares) agree to the observations to 
within a factor of three. The model results are again best fitted 
with two power laws, but the 13 CO observations do not extend 
to low enough N(H2) to confirm this. 

The observations show considerable scatter in the ratio be- 
tween A^CO) and A 7 ( 13 CO), with some sources deviating by a 
factor of two either way fro m the local elemental 12 C- 13 C ra- 
tio of 69 + 6 dWilsonll 19991) . It has long been assumed that the 
enhanced ratios are due to isotope-selective photodissociation 
and that the reduced ratios are due to isotope-exchange reac- 
tions. However, no answer has been found so far to the question 
of why one process do minates in some sources, and the other 
process in other sources JLanger et alj|1980t iMcCutcheon et af 
1 19801; ISheffer et al 



1992; Lambert et al. 1994; Federman et al 



20031; ISonnentrucker et al.M2007t ISheffer et alj|2007[) . Chemical 
models consistently show that isotope exchange is more effi- 
cient than selective photodissociation for observed cloud tem- 
peratures of up to ~60 K. The models can easily reproduce 
A 7 ( 12 CO)/A 7 ( 13 CO) ratios of less than 69, but they leave the ratios 
of more than 69 unexplained. 

Our grid of models suggests that the answer lies in supra- 
thermal chemistry. The right panel of Fig. [8] shows a plot of 
A 7 ( 12 CO)/A 7 ( 13 CO) versus A 7 (H 2 ) for the observations (grey) and 
for all depth steps at Ay < 5 from our models (black plus signs). 
Also shown are the 100-point means from the model, offset by a 
factor of three for clarity. The N( 12 CO) /A^ 13 CO) ratio initially 
increases from 69 to a mean of 151 at NQI2) = 4 x 10 20 cm" 2 , 
then decreases to a mean of 51 at 1 x 10 21 cm" 2 , and gradually 
increases again for larger depths. The turnover from the initial 
rise to the rapid drop is due to the transition from suprathermal 
to normal chemistry. 

At an Alfven speed of va = 3.3 km s _1 , the effective tem- 
perature for Reaction (0 is about 4000 K, so the forward and 
backward reactions are equally fast. Photodissociation is there- 
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photodissociation ceases to be the dominant destruction mechanism of I2 CO; the dotted line is the median. Likewise, the horizontal 
band indicates the range over which 12 CO becomes self- shielding. Right: ITJO-rjoint means of our model results (black squares), 
fitted with two straight lines . Also sh own are the fits to the observations bv lShefferet ail (2008. grey lines) and the results from the 
translucent cloud models o f|yDB88j (dashed lines) for* = 0.5, 1 and 10 {left to right). 
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fore the only active fractionation process and it enhances the 
amount of 12 CO relative to I3 CO. As we move beyond the depth 
to which the Alfven waves reach (N(H2) — 4 x 10 20 cirT 2 in 
our model), suprathermal chemistry is brought to a halt. The for- 
ward channel of Reaction (0 becomes faster than the backward 
channel and isotope exchange now reverses the fractionation. 
When we reach a depth of JV(H 2 ) = 1-2 x 10 21 cirT 2 , the 13 C + 
abundance drops too low to sustain further fractionation and the 



Af( 12 CO)/Af( 13 CO) ratio gradually returns to the elemental ratio 
of 69. 

The observations show a downward trend in the ratio be- 
tween JV( 12 CO) and JV( 13 CO) in the N(U 2 ) range from 2 x 10 20 to 
2 x 10 21 cirT 2 . This supports a rapid switch from selective pho- 
todissociation to isotope-exchange reactions being the dominant 
fractionation mechanism. A firm test of our model predictions 
requires observational data at both smaller and larger Af(H 2 ) than 
are currently available. 
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6.1.4. Other isotopologues 

Moving to the next two isotopologues, C n O and C ls O, the 
number of observed column densities goes down to three. 
lLambert et all dl994l) determined an 7V( 12 CO)//V(C 18 0) ratio of 
1550 ± 440 for the £ Oph diffuse cloud, some thre e times 
larger than the elemental 16 0- 18 ratio of 557 + 30 dWilsonl 
1999 ). An even hig her ratio of 3000 ± 600 was recorded by 
Shef fer et al.l d2002l) for the X Per translucent cloud, along with 
an /V( 12 CO)//V(C n O) ratio of 8700 ± 3600. The elemental I6 0- 
18 ratio is about four times smaller: 2000 ± 200. 

The H 2 column density towards X Per (HP 24534) is 8.5 x 
10 20 cm 2 (ISonnentrucker et all 120071; ISheffer et al.ll2007l) . For 
this value, our grid of models gives 7V( 12 CO)//V(C 18 0) = 680- 
3880 (median 1080) and 7V( 12 CO)/W(C 17 0) = 2520-10600 
(median 3660), consistent with the observations. Out of the three 
parameters that we varied, x affects the ratios most strongly. 
Both ratios decrease by a factor of two when x is increased 
from 1 to 10. The onset of self-shielding of 12 CO, and with 
it the onset of isotopic fractionation, occurs deeper into the 
cloud for a stronger UV field. Therefore, 7V( 12 CO)//V(C 18 0) and 
7V( 12 CO)/7V(C 17 0) are smaller at any given depth. 

In order to further explore the behaviour of C n O and C ls O, 
we take a more detailed look at our results for «h = 300 cm" 3 , 
Tgas = 50 K and^ = 1 (model C2ff from Table [T2l. This also al- 
lows for a comparison with 13 C n O and 13 C 18 Q. Both of these 
have been detecte d in interstellar clouds dLanger et al.l [1980; 
iBensch et al.ll200lT) . but column densities have not yet been de- 
rived. 

The photodissociation rates of CO and H2 for this combina- 
tion of parameters are plotted in Fig. [9] as a function of depth 
(z) into the cloud. The isotope-selective nature of the photodis- 
sociation is clearly visible. At a depth of 4.0 pc (Ay = 2.3 mag, 
N(U 2 ) = 1.7 x 10 21 cm" 2 ), 12 CO photodissociates about 7 times 
slower than 13 CO, 15-18 times slower than C 17 and C ls O, and 
23-26 times slower than 13 C n O and 13 C ls O. At depths between 
0.4 and 3.4 pc (Ay = 0.3-2.0 mag), the photodissociation of 
C ls O proceeds faster than that of C n O. The lines of the latter lie 
closer to those of 12 CO and are therefore more strongly shielded, 
as was the case for the diffuse £ Oph cloud (Table |4j. As the CO 
column grows larger in our model cloud, C ls O becomes self- 
shielding and its rate drops below that of C 17 around 3.4 pc. 

The extent to which each isotopologue is fractionated can 
easily be seen from the cumulative column density ratios when 
they are normalised to 12 CO and the isotope ratios: 



N z ( x C y O) [ 12 C] [ 16 0] 



N z ( 12 CO) [*C] pO] ' 
with [X] the elemental abundance of isotope X and 



(9) 



N z ( x aO)(z)= n( A C v O)dz' (10) 
Jo 

the cumulative column density. A plot of < R( 13 CO) as function 
of depth (Fig. ITOb shows the same trends as the model data in 
Fig. [8] the amount of 13 CO is reduced (ft < 1) up to a certain 
depth and is enhanced (ft > 1) further in by isotope-exchange 
reactions. We only see % < 1 for C n O and C ls O because 17 + 
and ls O + never become abundant enough to convert C lfi O in a 
reaction similar to Reaction d7). Photodissociation ceases to be 
the dominant destruction mechanism for C I7 and C 18 around 
Z = 4 pc (Fig.|9]l, so it can no longer cause substantial fractiona- 
tion and % gradually returns to unity. 

The remaining two isotopologues, I3 C 17 and 13 C ls O, expe- 
rience the I3 C and 17 0/ I8 fractionation simultaneously. Isotope 
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Fig. 9. Photodissociation rates of H2 and the six CO isotopo- 
logues as a function of depth into translucent cloud model C2a 
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the bottom axis indicate the depth at which He + takes over from 
UV radiation as the dominant destroyer for each isotopologue. 
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normalised to 12 CO and the elemental isotopic abundances (Eq. 
(O) as a function of depth into translucent cloud model C2a. 



exchange enhances their abundances relative to C 17 and C 18 
in the same manner and at the same depths as it enhances the 
abundance of 13 CO relative to 12 CO. However, isotope-selective 
photodissociation is a stronger effect for the two heaviest iso- 
topologues, so % remains less than unity. This result holds across 
our entire grid of translucent cloud models. 
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6.2. Photon-dominated regions 

Photon-dominated regions (PDRs) form another type of objects 
where photodissociation is a key process in controlling the CO 
abundance. Their densities are higher («h up to ~10 6 cm 4 ) than 
those of diffuse and translucent cl ouds and they are exposed t o 
stronger UV fields Of up to ~10 6 ; iHollenbach & Tielenslll997h . 
Observational efforts have mostly gone into mapping the struc- 
ture of PDRs, so CO column densities are tabulated only for a 
small number of sources. For most of these the column density 
of neutral carbon has also been determined. There are no direct 
measurements of H2 column densities; instead, 7V(H2) is usually 
obtained from jV( 12 C O) using a typic al abundance ratio between 
the two species (e.g., Frerking et al.ll 1982b . 

We present here the results from our combined photodisso- 
ciation and chemistry model for a grid of physical parameters 
appropriate for PDRs: r gas = 50, 80 and 100 K, « H = 10 3 , 10 4 , 



10 5 and 10 b 



and x — 1 , 10 and 10 . The excitation 



temperatures of CO and H2 are set equal to r gas . Where nec- 
essary, we expand the number of rotational levels of CO beyond 
the default limit of J"-l, and likewise for H 2 (see Sect.l4~4l>. UV 
pumping increases the population in the v"=0, J">3 and v">0 
levels of H2. However, this does not affect the CO photodissoci- 
ation rates at these temperatures, so the level population of H2 is 
simply taken as fully thermal. The models are run to an Ay of 30 
mag and are otherwise unchanged from the previous section. 

Figure [TT| shows the relationship between the column den- 
sities of 12 CO and 12 C for a large series of depths steps from 
our grid. The data for «h = 10 4 and 10 6 cm 4 are omitted for 
reasons of clarity. Overplotted are column densiti es determined 
towards a number of posit i ons in different P DRs (Beutheret al. 
2000; Kamega reTai.ll2003l;lRramer et alj2008HSun et alj|2 008). 
Figure 5 a of lMookeriea et al] d2006l) contains additional data for 
several dozen positions in the PDRs in Cepheus B. In none of 
these cases was the 12 CO column density measured directly; 
instead, it was derived from the 13 CO or C ls O column den- 
sity using standard abundance ratios. The spread in both the 
model predictions and the observational data is large. The col- 
umn of 12 C for each individual model is nearly constant for 
7V( 12 CO) > 3 x 10 17 cirT 2 , as shown by the solid line connect- 



ing the various depth steps for nu 



gas 



50 K 



and^ = 10 . The observations seem to correspond more to the 
lower-density end of the model grid (the plus signs), but the sam- 
ple is too small and is too scattered to be conclusive. 

The onset of shielding is seen to occur around A^( 12 CO) = 
10 15 cirT 2 in the models. From that point onwards, the low- 
density models gain more atomic carbon than do the high- 
density models. For a given column of 12 CO, the photodisso- 
ciation rates (or, in other words, the 12 C formation rates) are 
roughly independent of density, but the destruction of 12 C oc- 
curs faster in higher densities. The same argument explained the 
density dependence of the N( l2 CO)/N(H2) ratio in the translu- 
cent clouds. 

The two isotopologues for which col umn densities have been 
determined directly are C 1 8 Q and 13 CO dLorenll989tlGerin et alj 
19981: iPlume et all 119991: iBeuther et alj 120001; ISchneider et al] 
2003[ iMookeriea et alJ I2OO6I: iKramer etal J 120081; ISunetalJ 
2008). Again, both the observations and the models show a large 



spread in column densities (Fig. [T2l . It seems the observations 
do not trace the high-extinction (Ay > 10 mag), high-density 
(«h > 10 5 cm" 3 ) material shown in the upper left corner of each 
panel, but a more detailed study is required to draw any firm 
conclusions on this point. N( l2 CO), 7V(C ls O) and N( n CO) are 
mutually well correlated in both the observations and the mod- 



els, but significant deviations can easily occur at specific depths 
due to the different photodissociation rates. When using standard 
abundance ratios to derive the column density of 12 CO from that 
of another isotopologue, as was done for all cases in Fig.Ql] the 
result will generally be accurate to at best a factor of two. 

6.3. Circumstellar disks 
6.3.1. Model setup 

Low-mass stars like our own Sun are formed through the grav- 
itational collapse of a cold molecular cloud. As the young star 
grows at the centre, it gathers part of the remaining cloud ma- 
terial into a circumstellar disk. It is inside this disk that planets 
are formed, so its chemical composition is of great interest. The 
physical structure of the disk may be simplified as consisting of 
a cold, dense region near the midplane, covered by a warmer 
region of lower density. This surface layer intercepts the star's 
radiation and shows many similarities t o common PDRs . 

We adopt the standard disk model of lD'Alessio et afljl999), 
whose chemistry has been studied extensively ( Aikawa et alj 
l2002Hvan Zadelhoffet alJl2003b I.Tonkheid et alJl2004 . The star, 
a T Tauri type, has a mass of 0.5 M , a radius of 2 R Q and an 
effective temperature of 4000 K. It is surrounded by a disk of 
mass 0.07 M Q and outer radius 400 AU. The disk is accreting 
onto the star at a constant rate of 10~ 8 M yr~' and its viscosity 
is characterised by a = 0.01 . We focus on a vertical slice through 
the disk at a radius of 105 AU, where the surface is located at a 
height of z = 120 AU. This slice is irradiated from the top by a 
radiation field with an effective UV intensi ty of 516 times that 
of the interstellar radiation field. Fo llowing Ivan Zadelhoff et afl, 
we adopt two spectral shapes: the iDraine field and a 4000 K 
blackbody spectrum. The latter is very weak in the wavelength 
range where CO is photodissociated. We adopt again the chem- 
ical network from Sect. 16.1.11 and we set the gas te mperature 
equal to the dust temperature from the iD'Alessio et alj model. 
Given the high densities (2 X 10 4 cm~ 3 at the surface and more 
than 10 9 cirT 3 at the midplane), we also use r gas for r ex (CO) and 
TexiHi)- As in the PDR models, UV pumping of H2 does not af- 
fect the CO photodissociation rates, so the H2 level populations 
are taken as fully thermal. 

In reality, disks are not irradiated by the iDraind field or 
a pure blackbody. For example, T Tauri stars typically show 
UV emission in excess of these simple model spectra. Many of 
them also emit X-rays f rom accretion hotspots. As shown by 
Ivan Zadelhoff et alj d2003l) . the photochemistry of the disk can 
be sensitive to the details of the radiation field. Moreover, the 
excess emission can be concentrated in a few specific lines such 
as Hi Lyman-ff. CO cannot be dissociated by Ly-a, but it can 
be by the Hi Ly-y line at 972.54 A through overlap with the W0 
band at 972.70 A (No. 25 from TableQ}- A detailed treatment of 
these effects is beyond the scope of the current work. 



6.3.2. Isotopologue ratios 

The depth depend ence of the photodissociation rates in the disk, 
illuminated by the lDraing field, is qualitatively the same as in the 
cloud and PDR models from Sects. 16. II and [6721 Several quan- 
titative differences arise because of the different UV intensity 
and the higher densities. The profiles of the normalised column 
density ratios, % (Eq. (|9j), do not show any enhancement in 
13 CO. Moving down from the surface, suprathermal chemistry 
initially prevents the formation of 13 CO from 12 CO and 13 C + . 
The Alfven waves penetrate to 7V(H2) 
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Fig. 11. Column density of 12 CO versus that of 12 C on two different scales. Black: results from our grid of PDR models with 
«h = 10 3 cirT 3 (plus signs) and 10 5 cuT 3 (diamonds). The line traces the re s ults for differen t depth steps at »h = 10 s c m -3 , 
r gas = 50 K an d x = 10 4 - Grey squares: observations of iBeuther eTall (l2000h . iKamegai etall (l2003h . iKramer et aT] (120081) and 
Sun etall (120081) . 
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Fig. 12. Column densities of C ls O (left) and I3 CO (right) versus that of I2 C. Symbols are as in Fig.Qj] See text for references to 
the observations. 



translates here to a height of about 50 AU. The absence of supra- 
thermal chemistry below that point allows < R( 13 CO) to increase 
again. However, the 13 C + abundance decreases at the same time, 
and 7?( 13 CO) never gets above unity. This is consistent with the 
/V( 12 CO)//V( 13 CO) ratios of more than 69 determined for the 
disks around HL Tau AB Aur and VV CrA dBrittain et alj2005t 
ISchrever et aDl2008t ISmith et al.ll2009h . 

In terms of loca l abundances, our results s how a trend oppo- 
site to the model by Woods & Willacv (2009). Moving up verti- 
cally from the midplane at a given radius, their n( 12 CO)/n( 13 CO) 
abundance ratio initially remains nearly constant, then drops 
to a lower value, and finally shows a shallow rise for the re- 
mainder of the column. We also find a constant ratio at first, 
but then we get a rise (equivalent to 7?( 13 CO) < 1) and a 



gradual decrease back to the elemental ratio of 69. The dif- 
fer ence is due to th e absence of suprathermal chemistry in 
the|Woods & Willacv model. Suprathermal chemistry effectively 
shuts down the isotope-exchange reactions, leaving photodisso- 
ciation (which increases «( 12 CO)/«( 13 CO)) as the only fraction- 
ation mechanism. Additional observations of 12 CO and 13 CO in 
circumstellar disks are required to ascertain whether isotope ex- 
change or selective photodissociation is the key process in con- 
trolling the ratio of the abundances of these two species. 

For C n O and C ls O, H is ~1 from the disk surface down to 
70 AU. It then gradually drops to ~0.1 due to their photodissoci- 
ation rates being higher than that of 12 CO. Shielding of C n O and 
C 18 sets in around 32 AU, so that photodissociation no longer 
causes isotope fractionation and 7? increases again. However, 



18 



R. Visser et al.: The photodissociation and chemistry of interstellar CO isotopologues 



Am (mag) 



10 2 10 1 10° 10 




10 2 10 1 10° 10" 



0.5 - 



0.0 



a gr = 1 jj,m 



O ' 1 13 C0 



100 120 
z (AU) 




17 0, C 18 0, 13 C 17 0, 13 C 18 



80 



100 



120 



Fig. 13. Cumulative column densities of the CO isotopologues normalis ed to 12 CO and the elemental isotopic abundances (Eq. (0) 
as a function of height at a radius of 105 AU in our disk model with the Draine field and dust grains of 0. 1 yum (left) or 1 //m (right). 
CO freezes out below 24 AU. The diamond symbols indicate where H3 takes over from UV photons as the main destroyer of each 
isotopologue. 



it does not reach a value of unity. Below 24 AU, the temper- 
ature is low enough that most of the CO freezes out onto the 
dust grains. The tiny amount remaining in the gas phase does 
not add to the column density, so % is effectively frozen at 0.75 
and 0.73 (Fig. [TBI. The corresponding column density ratios are 
AT( 12 CO),W(C 17 0) = 2700 and jV( 12 CO),W(C 18 0) = 760. Both 
values are consistent with the degree of fracti onation determine d 
for the disk around the T Tauri star VV CrA dSmifh et al.ll2009h . 

6.3.3. Model variations 

If we illuminate the disk by a 4000 K blackbody (spectrum C 
from van Zadelhoff et al. 2003), it does not receive enough pho- 
tons in the 912-1118 A range to cause substantial photodisso- 
ciation of CO. Instead, the non-isotope-selective reactions with 
He + and are the main destruction mechanisms even at the 
top of our vertical slice at R — 105 AU. The normalised column 
density ratios therefore remain unity at all heights. Adding the 
interstellar UV field at% = 1 to the 4000 K blackbody spectrum 
allows CO to be photodissociated again in the surface layers. 
The isotope fractionation returns, but only partially, because the 
UV field is still weaker than in the original case. Freeze-out now 
locks <R(C n O) and <R(C 18 0) at 0.86 and 0.85 instead of -0.75. 
Clearly, some UV component is required in the irradiating spec- 
trum to get isotope fractionation through the photodissociation 
of CO, but it does not have to be a very strong component: the 
interstellar UV field already has a significant effect in this model. 

Another parameter to be varied is the grain size. The 
disk model so far contained small dust grains, with an av- 
erage size of 0.1 jum. Grain growth is an inevitable step to- 
wards planet formation and is indeed known to occur in disks 
(Ivan Boekel et al.l 120031: iPrzygodda et all 12003b lLommen etai] 
2007). Larger grains cause less extinction in the UV dShen et all 
2004), so CO can be photodissociated deeper into the disk. Grain 
growth also affects the H2 formation rate: the rate is inversely 
proportional to app roximately the square of the grain radius 
dJonkheicT et al. 2006). The right panel of Fig.[T3"lshows the nor- 
malised vertical column density ratios at a radius of 105 AU for 
an average grain size of 1 yum. The onset of fractionation occurs 
at a lower altitude than in the 0. 1 jum case: 55 AU instead of 75 



AU. The slower formation rate of H2 has a direct influence on 
the formation of CO through the pathway starting with Reaction 
©. There is less H2 and CO in the disk's surface layers when 
we increase the grain size, so we have to go down deeper before 
12 CO becomes self-shielding and isotope fractionation sets in. 
At the same time, the weaker continuum extinction allows the 
UV radiation to penetrate deeper, so the isotope-selective pho- 
todissociation goes on to lower altitudes. The spike in "R( 13 CO) 
to 1 .4 at 30 AU is due to a narrow zone where the temperature 
is relatively low (35 K), the Alfven waves do not reach, and the 
13 C + abundance is high enough to convert 12 CO into 13 CO. The 
model with 0.1 //m grains also has a narrow zone with a high 
13 C + abundance and no suprathermal chemistry, but the kinetic 
temperature in this zone is 50 K. There is still some conversion 
of 12 CO into 13 CO at that temperature, but not enough to show 
as a spike in the normalised column density ratio. 

The bigger dust grains also affect the fractionation of C 17 
and C 18 0. The point at which photodissociation ceases to be the 
main destruction mechanism now nearly coincides with the point 
where CO freezes out (Fig. [T3l . Both species have % « 0.15 
when this happens, so the ratios are locked at that value. In the 
small-dust model, the ratios were locked at ~0.75. Observations 
of face-on disks in principle probe CO all the way to the mid- 
plane, so our model predicts that grain growth from 0.1 to 1 //m 
can increase the observed fractionation of C n O and C 18 by a 
factor of five. The fractionation in the local abundances, which 
is relevant for incorporating material into comets, increases by 
the same amount from 0.1 to 1 //m grains. 

All the results so far in this section are for a T Tauri disk. 
Herbig Ae/Be stars are more massive and more luminous and 
tend to have warmer disks, where CO may not be frozen out at 
all. That would prevent % from being locked as it is in Fig. |T3l 
instead, it would rapidly return to unity for all isotopologues. 
Observations probing the midplane of such a disk would not find 
any fractionation in n O and 18 relative to 16 0. 

6.3.4. Implications for meteorites 

Regardless of the grain size and the absence or presence of 
freeze-out, the model shows n O and 18 to be nearly equally 
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fractionated: % has very similar values for C 17 and C ls O. The 
equal photodissociation rates implied by this result are partly 
due to our choice of molecular parameters. Spectroscopic data 
on C I7 are still very scarce, so we assume many of its oscilla- 
tor strengths and predissociation probabilities to have the same 
value as for C ls O (Sects. |231 and l2~4l). Still, two important dif- 
ferences remain. First, the elemental abundance of 18 is higher 
than that of 17 0, resulting in some self-shielding for C I8 (Fig. 
[TJ- Second, the shifts of the C 17 lines with respect to 12 CO 
are generally somewhat smaller than those of the C 18 lines, so 
C I7 is more strongly shielded by I2 CO. Both effects are minor, 
however, and the photodissociation rates end up being nearly the 

same. 

According to the experiments of Chakrabo rtv et al.l (120081) . 
the photodissociation rate due to the EQ and El bands may be 
up to -40% higher for C n O than for C 18 0. If we mimic this 
by artificially increasing the E0 and El oscillator strengths of 
C 17 by 40%, the overall change in 7?(C n O) is less than 10% at 
any depth. When the two oscillator strengths are doubled from 
their standard value, 7?(C 17 0) changes by at most 20%. In both 
cases, the change in 7?(C 17 0) is largest at the top of the disk and 
gradually decreases for lower altitudes; it is less than 1% below 
50 AU. 

Observational evidence for equa l fractionation of 1 7 Q and 
18 in disks was recently obtained by Smit h et all d2009l) for the 
T Tauri star VV CrA. Their analysis points to isotope-selective 
photodissociation of CO as the most likely explanation, although 
they cautioned that more observations are needed to confirm that 
conclusion. Meanwhile, their observations are consistent with 
our model results. 

A long-standing puzzle in our solar system involves the 
equal fractionat ion of 17 Q and 18 Q in the most refractory phases 
of meteorites dClavton et al.l Il973l) . Current evidence points 
at isotope-selective phot odissocation o f CO in the solar neb- 
ula playing a key role dClavtonl l2002t iLvons & Younjd l2005; 
iLee et al.l2008l) . It is generally accepted that the resulting atomic 
oxygen eventually makes it into the refractory material to pro- 
duce the observed isotope ratios, but it remains unknown how 
this actually happens. Any additional mass-independent frac- 
tionation during that process must be limited, because pho- 
todissociation of C O already rep roduces the observed isotope 
abundance ratios dYoungl 120071) . One important requirement 
here is to have C n O and C ls O photodissociate at the same 
rate, an assumption that had gone untested up to now. Our 
model confirms they do so in many different environments — 
including diffuse and translucent clouds, PDRs, and circumstel- 
lar disks — even when allow ing for the isotope effects suggested 
bv lChakrabortv et all d2008l) . 

Equal (or nearly equal) 17 and ls O fractionation also 
occurs for other oxygen-containing molecules in the model, 
such as water. The equal fractionation is a clear sign of CO 
photodissociation playing a key role in the fractionation pro- 
cess. Other fractionation processes, such as isotope-exchange 
reactions, diffusion, freeze-out and evaporation, depend on the 
vibrational energies of the isotopologues and are therefore 
mass-dependent. Specifically, such processes would lead to 
17 Q b eing 0.52 times as fractionated as 18 (Matsuhisa et alj 
fl978h . If mass-independent fractionation processes other than 
isotope-selective photodissociation play a role in circumstel- 
lar disks, they have not been discovered so fa r ( Cravton| [2002; 
Yurimoto & Kuramotoll2004l: ILvons & Youndl2005l) . Recently. 
Kimura et al. (2007) did report on a set of chemical experiments 
yielding mass-independently fractionated silicates, but they were 
unable to identify the underlying reaction mechanism. Hence, 



it remains unkn own if this would be relevant for disks. Self- 
shielding in O2 dThiemens & Heidenreich|[l983l) is unimportant 
because O2 never becomes abundant enough to become a signif- 
icant source of atomic O. 



7. Conclusions 

This paper presents an update d model, based on the method of 
van Dishoeck & Black ( 1988), for the photodissociation of car- 
bon monoxide (CO) and its heavier isotopologues. It contains 
recent spectroscopic data from the literature and produces a pho- 
todissociation rate of 2.6 x 10~ 10 s in the interstellar medium, 
30% higher than currently adopted values. Our model is the first 
to include C n O and 13 C n O and we apply it to a broaderrange of 
model parameters than has been done before. The main results 
are as follows: 

• Self-shielding is very important for 12 CO and somewhat 
so for C ls O and 13 CO. The rare isotopologues 13 C n O 
and 13 C ls O have the highest photodissociation rate at most 
depths into a cloud or other object. The rates of C 17 and 
C ls O are very similar at all depths, but details in the line 
overlap and self-shielding can cause mutual differences of 
up to 30% (Sects.|331and |6X4l 

• When coupled to a chemical network, the model reproduces 
column densities observed towards diffuse clouds and PDRs. 
It shows that the large spread in observed 7V( 12 CO)/7V( 13 CO) 
ratios may be due to a combination of isotope-selective pho- 
todissociation and suprathermal chemistry (Sect. l6"T31 ). 

• Photodissociation of CO is more strongly isotope-selective 
in cold gas than in warm gas (Sect.[4]i. 

• The results from the full computation are well approximated 
by a grid of pre-computed shielding functions, intended for 
easy use in various other models. Shielding functions are 
provided for a range of astrophysical parameters (Sect. [5]). 

• Grain growth in circumstellar disks increases the vertical 
range in which CO can be photodissociated. If photodisso- 
ciation is still important at the point where CO freezes out 
onto the grains, the observable gas-phase column density ra- 
tios 7V( 12 CO)/7V(C 17 0) and 7V( 12 CO)/7V(C 18 0) may become 
an order of magnitude larger than the initial isotopic abun- 
dance ratios (Sect. 16331 ). 

• Without a far-UV component in the irradiating spectrum, 
isotope fractionation through the photodissociation of CO is 
not possible. The interstellar radiation field can already cause 
substantial fractionation in disks (Sect. [633]). 

• The isotope-selective nature of the CO photodissociation 
results in mass-independent fractionation of n O and 18 0. 
Column density ratios computed for a circumstellar disk 
agree well with recent observations. Our model supports the 
hypothesis that the photodissociation of CO is responsible 
for the anomalous 17 and 18 abundances in meteorites 
(Sect. |6331). 
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Table 1. New molecular data for 12 CO. a 
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Re ferences : ( 1 ) ^ and v from lUbachs et alj (Il994l): (2) A ,„, from lUbachs et alj (tl994j ); (3) B[, and D' fromlUbachs et all dl994l): (4) An an d v n fromlvDB88l: (5 ) B[, and D[, from 



lyDB88l: (6) f v0 and A tot from lEidelsberg et alj d2006l): ( 7) A , v , B' w an d D[, fromlCacciani et alj d2002l) : (8) A t „ t fromlDrabbels et all dl993l) : (9) f v0 fro mlEidelsberg etal.1 
d2004 : (10) M h v , B' v and D'i from lUbachs et alj d2000): (11) f v , from Eidelsberg^eLaU d2006j); (12) A tot from lUbachs et all d2000l) ; (13) cj' e , oj e x' e from lKepaldl988l) : (14) 



<k>, vo, Atot, B[, and £>.'. from ICacciani et al.1 (1200 ll): (15) /v^from iFederman et all (l200ll) : (16) tj' e , oj e x' e from lTilford & Vandershcel (Il968l) : (17) A Q , v , B' v and D' v from 



ICacciani & Ubachsl (12004 : (181 A tot from lCacciani et al. (1998) 
Notes: 

a /lo, vo, ID and v' from Eidelsbe rg et alJ ([1992), / v <o, A to t, B' v and D' v from ER90, unless noted otherwise. All rj are recomputed from the f v >o and A to t- The notation a(b) in this 

and following tables m eans a X 10 b . Many values are rounded off from higher-precision values in the references; 

b The numbering follows lvDB88l Their bands 1 and 2 are split in to four and two components. The corresponding ER90 indices are also given. 
c a>' e and a> e x' e from the CO ground state (iGuelachvili et alj|1983l) . ca e y' e and u> e z' e (not listed) are included in the model. 
d cj' s and w e x; from the CO + X 2 E + sta te faaridass et alj|2000^ 
e B' v and D' v from the CO + A 2 II state dHaridass et alJl2000h. 

f a>' e and <jj e jCg from the CO + A 2 FL state (Haridass et al. 2000). Lo e y' e (not listed) is included in the model. 

g Vibrational constants derived from th e different vp in one o f six vibrational series: bands 30-27, 28-21, 25-20-16-8, 24-19, 22-18 or 15-5. 
h B' v and D' w from the CO + X 2 Z + state dHaridass et al.ll2000h . 

1 A tot depends on parity and/or rotational level (see Ta ble [3]l. A to t and rj are listed here for J'-O and / parity, 
j B'„ and D'., computed from the C 18 Q values of lER90l 

k B' v and D' v depend on parity (see Tabled; values are listed here for / parity. 
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Table 6. As Tableg] for fo(CO) = 0.3 km s"\ 7/ ex (CO) = 50 K and N( 12 CO)/7V( 13 CO) = 69. a 



logMH 2 ) log 12 CO (cm- 2 ) 

(cm' 2 ) 13 14 15 16 17 18 19 

12 CO: unattenuated rate k 0J = 2.590 x KT 10 s' 1 

L000 9.405(-l) 7.046(-l) 4.015(-1) 9.964(-2) 1.567(-2) 3.162(-3) 4.839(-4) 

19 7.546(-l) 6.979(-l) 4.817(-1) 2.577(-l) 6.505(-2) 1.135(-2) 2.369(-3) 3.924(-4) 

20 5.752(-l) 5.228(-l) 3.279(-l) 1.559(-1) 3.559(-2) 6.443(-3) 1.526(-3) 2.751(-4) 

21 2.493(-l) 2.196(-1) 1.135(-1) 4.062(-2) 7.864(-3) 1.516(-3) 4.448(-4) 9.367(-5) 

22 1.550(-3) 1.370(-3) 6.801(-4) 2.127(-4) 5.051(-5) 1.198(-5) 6.553(-6) 3.937(-6) 

23 8.492(-8) 8.492(-8) 8.492(-8) 8.492(-8) 8.492(-8) 8.492(-8) 8.488(-8) 8.453(-8) 
C' 7 Q: unattenuated rate k 0J = 2.607 x 10~'° s~' 

L000 9.756(-l) 8.826(-l) 7.413(-1) 4.507(-l) 1.508(-1) 2.533(-2) 1.684(-3) 

19 7.266(-l) 7.036(-l) 6.209(-l) 5.209(-l) 3.212(-1) 1.074(-1) 1.832(-2) 1.436(-3) 

20 5.418(-1) 5.206(-l) 4.469(-l) 3.695(-l) 2.222(-l) 7.000(-2) 1.197(-2) 9.425(-4) 

21 2.208(-l) 2.089(-l) 1.687(-1) 1.352(-1) 8.073(-2) 2.082(-2) 3.360(-3) 2.492(-4) 

22 1.401(-3) 1.352(-3) 1.190(-3) 1.043(-3) 7.022(-4) 1.687(-4) 2.525(-5) 4.21 1(-6) 

23 8.509(-8) 8.509(-8) 8.509(-8) 8.509(-8) 8.509(-8) 8.509(-8) 8.505(-8) 8.470(-8) 
C' 8 Q: unattenuated rate k 0J = 2.396 x 10~'° s~' 

L000 9.822(-l) 9.1630-1) 8.067(-l) 5.498(-l) 2.188(-1) 3.412(-2) 1.992(-3) 

19 8.007(-l) 7.833(-l) 7.206(-l) 6.272(-l) 4.226(-l) 1.597(-1) 2.469(-2) 1.630(-3) 

20 5.848(-l) 5.688(-l) 5.123(-1) 4.363(-l) 2.816(-1) 9.826(-2) 1.653(-2) 1.158(-3) 

21 2.277(-l) 2.187(-1) 1.881(-1) 1.546(-1) 9.811(-2) 3.114(-2) 5.506(-3) 3.655(-4) 

22 1.411(-3) 1.375(-3) 1.256(-3) 1.126(-3) 8.200(-4) 2.912(-4) 5.310(-5) 5.241(-6) 

23 8.937(-8) 8.937(-8) 8.937(-8) 8.937(-8) 8.937(-8) 8.936(-8) 8.933(-8) 8.896(-8) 
13 CO: unattenuated rate k 0J = 2.597 x 10~'° s~' 

L000 9.765(-l) 8.965(-l) 7.701(-1) 4.459(-l) 1.415(-1) 1.748(-2) 8.346(-4) 

19 7.780(-l) 7.550(-l) 6.789(-l) 5.753(-l) 3.262(-l) 1.032(-1) 1.183(-2) 6.569(-4) 

20 5.523(-l) 5.309(-l) 4.617(-1) 3.793(-l) 1.958(-1) 5.888(-2) 7.515(-3) 4.382(-4) 

21 2.100(-1) 1.979(-1) 1.601(-1) 1.244(-1) 5.373(-2) 1.573(-2) 2.535(-3) 1.361(-4) 

22 1.318(-3) 1.268(-3) 1.107(-3) 9.017(-4) 3.658(-4) 1.114(-4) 2.414(-5) 2.608(-6) 

23 4.511(-8) 4.511(-8) 4.511(-8) 4.511(-8) 4.511(-8) 4.511(-8) 4.509(-8) 4.490(-8) 
13 C 17 Q: unattenuated rate fc 0- , = 2.605 x lO' 10 s"' 

1.000 9.853(-l) 9.344(-l) 8.4530-1) 5.9780-1) 3.075(-l) 1.180(-1) 3.205(-2) 

19 8.097(-l) 7.954(-l) 7.474(-l) 6.7550-1) 4.9470-1) 2.7110-1) 1.065(-1) 2.898(-2) 

20 5.925(-l) 5.7920-1) 5.363(-l) 4.8150-1) 3.5760-1) 2.062(-l) 8.289(-2) 1.967(-2) 

21 2.389(-l) 2.315(-1) 2.083(-l) 1.848(-1) 1.410(-1) 8.627(-2) 3.1850-2) 3.489(-3) 

22 1.937(-3) 1.907(-3) 1.8120-3) 1.71 1(-3) 1.4740-3) 1.042(-3) 3.690(-4) 1.704(-5) 

23 4.523Q-8) 4.523Q-8) 4.523Q-8) 4.523Q-8) 4.523Q-8) 4.523Q-8) 4.521Q-8) 4.502Q-8) 
13 C 18 Q: unattenuated rate k 0ii = 2.510 x 1Q-' s" 1 

L000 9.873(-l) 9.424(-l) 8.744(-i) 6.740(-l) 3.8850-1) 1.2680-1) 2.980(-2) 

19 7.9750-1) 7.852(-l) 7.424(-l) 6.8940-1) 5.4470-1) 3.236(-l) 1.036(-1) 2.722(-2) 

20 5.8690-1) 5.754(-l) 5.364(-l) 4.953(-l) 3.932(-l) 2.445(-l) 8.070(-2) 1.894(-2) 

21 2.4650-1) 2.401(-1) 2.1870-1) 2.01 1(-1) 1.6280-1) 1.079(-1) 3.267(-2) 4.406(-3) 

22 2.164(-3) 2.1380-3) 2.050(-3) 1.971(-3) 1.772(-3) 1.4210-3) 5.147(-4) 4.066(-5) 

23 4.454(-8) 4.454(-8) 4.454(-8) 4.454(-8) 4.454(-8) 4.454(-8) 4.452(-8) 4.434(-8) 



These shielding functions were computed for the same parameters as in Table|5] except r cx (CO) = 50 K and r cx (H2) = 50 1 5 K. Additional 
rotational lines of CO and IT were included as described in Sect. 14.41 
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Table 7. As Tableg] for b(CO) = 3.0 km s"\ r ex (CO) = 5 K and JV( 12 CO)/jV( 13 CO) = 69. a 



logMH 2 ) log 12 CO (cm- 2 ) 

(cm' 2 ) 13 14 15 16 17 18 19 

12 CO: unattenuated rate k 0J = 2.588 x IP' 10 s~' 

L000 9.639(-l) 7.443(-l) 3.079(-l) 5.733(-2) 1.321(-2) 3.322(-3) 5. 193 (-4) 

19 7.928(-l) 7.594(-l) 5.606(-l) 2.145(-1) 4.172(-2) 1.039(-2) 2.696(-3) 4.426(-4) 

20 7.037(-l) 6.743(-l) 4.984(-l) 1.861(-1) 3.410(-2) 8.272(-3) 2.212(-3) 3.716(-4) 

21 3.176(-1) 3.082(-l) 2.477(-l) 1.007(-1) 1.833(-2) 4.252(-3) 1.066(-3) 1.850(-4) 

22 1.083(-2) 1.067(-2) 9.413(-3) 4.528(-3) 1.054(-3) 2.353(-4) 6.789(-5) 1.378(-5) 

23 3.931(-7) 3.931(-7) 3.930(-7) 3.922(-7) 3.896(-7) 3.886(-7) 3.883(-7) 3.869(-7) 

C' 7 Q: unattenuated rate k 0J = 2.604 x 1(T' s~' 

L000 9.703(-l) 7.959(-l) 5.192(-1) 3.667(-l) 2.247(-l) 6.359(-2) 6.394(-3) 

19 8.049(-l) 7.764(-l) 6.108(-1) 3.734(-l) 2.733(-l) 1.692(-1) 4.882(-2) 5.903(-3) 

20 7.008(-l) 6.759(-l) 5.31 1(-1) 3.212(-1) 2.348(-l) 1.424(-1) 3.533(-2) 3.847(-3) 

21 3.125(-1) 3.056(-l) 2.638(-l) 1.885(-1) 1.474(-1) 8.575(-2) 1.807(-2) 1.839(-3) 

22 1.077(-2) 1.076(-2) 1.062(-2) 9.872(-3) 8.416(-3) 5.126(-3) 1.119(-3) 4.151(-5) 

23 3.950(-7) 3.950(-7) 3.950(-7) 3.950(-7) 3.949(-7) 3.945(-7) 3.919(-7) 3.864(-7) 

C' 8 Q: unattenuated rate k 0J = 2.394 x 10~'° s~' 

LOOO 9.767(-l) 8.397(-l) 6.18K-1) 4.521(-1) 1.999(-1) 3.383(-2) 1.883(-3) 

19 8.114(-1) 7.886(-l) 6.554(-l) 4.491(-1) 3.245(-l) 1.475(-1) 2.665(-2) 1.723(-3) 

20 7.177(-1) 6.979(-l) 5.816(-1) 3.981(-1) 2.851(-1) 1.284(-1) 2.244(-2) 1.476(-3) 

21 3.313(-1) 3.260(-l) 2.936(-l) 2.274(-l) 1.646(-1) 7.081(-2) 1.188(-2) 6.966(-4) 

22 1.130(-2) 1.129(-2) 1.113(-2) 1.016(-2) 7.643(-3) 3.131(-3) 6.377(-4) 3.195(-5) 

23 4.233(-7) 4.233(-7) 4.233(-7) 4.233(-7) 4.231(-7) 4.217(-7) 4.167(-7) 4.126(-7) 
13 CO: unattenuated rate k 0J = 2.595 x 10~'° s~' 

LOOO 9.750(-l) 8.292(-l) 6.083(-l) 4.222(-l) 1.249(-1) 1.250(-2) 6.838(-4) 

19 7.899(-l) 7.656(-l) 6.246(-l) 4.244(-l) 2.974(-l) 9.238(-2) 9.445(-3) 5.851(-4) 

20 6.933(-l) 6.720(-l) 5.484(-l) 3.694(-l) 2.561(-1) 7.843(-2) 7.799(-3) 4.678(-4) 

21 3.278(-l) 3.219(-1) 2.858(-l) 2.162(-1) 1.496(-1) 4.297(-2) 4.117(-3) 2.360(-4) 

22 1.108(-2) 1.106(-2) 1.089(-2) 9.813(-3) 6.942(-3) 2.057(-3) 2.618(-4) 1.306(-5) 

23 2.378(-7) 2.378(-7) 2.378(-7) 2.378(-7) 2.375(-7) 2.354(-7) 2.304(-7) 2.285(-7) 
13 C 17 Q: unattenuated rate fc 0- , = 2.604 x 10''° s' 1 

1.000 9.800(-l) 8.605(-l) 6.553(-l) 5.067(-l) 2.970(-l) 1.491(-1) 4.514(-2) 

19 8.134(-1) 7.939(-l) 6.781(-1) 4.891(-1) 3.826(-l) 2.401(-1) 1.228(-1) 3.233(-2) 

20 7.049(-l) 6.879(-l) 5.867(-l) 4.192(-1) 3.270(-l) 2.096(-l) 1.079(-1) 2.669(-2) 

21 3.302(-l) 3.256(-l) 2.970(-l) 2.382(-l) 1.931(-1) 1.292(-1) 6.607(-2) 1.136(-2) 

22 1.174(-2) 1.173(-2) 1.161(-2) 1.093(-2) 9.517(-3) 7.132(-3) 3.652(-3) 3.436(-4) 

23 2.312(-7) 2.312(-7) 2.312(-7) 2.312(-7) 2.31 1(-7) 2.31 1(-7) 2.308(-7) 2.290(-7) 
13 C 18 Q: unattenuated rate k 0ii = 2.509 x 10''° s" 1 

L000 9.826(-l) 8.767(-l) 6.711(-i) 5.145(-1) 3.329(-l) 1.592(-1) 4.547(-2) 

19 8.157(-1) 7.988(-l) 6.964(-l) 5.074(-l) 3.964(-l) 2.777(-l) 1.330(-1) 4.212(-2) 

20 7.096(-l) 6.950(-l) 6.057(-l) 4.394(-l) 3.440(-l) 2.477(-l) 1.197(-1) 3.796(-2) 

21 3.192(-1) 3.153(-1) 2.911(-1) 2.365(-l) 1.943(-1) 1.424(-1) 6.786(-2) 1.888(-2) 

22 1.077(-2) 1.076(-2) 1.065(-2) 9.980(-3) 8.532(-3) 6.289(-3) 2.998(-3) 6.921(-4) 

23 2.410(-7) 2.410(-7) 2.410(-7) 2.410(-7) 2.410(-7) 2.409(-7) 2.407(-7) 2.395(-7) 



These shielding functions were computed for the same parameters as in Table \5\ except b(CO) = 3.0 km s b(H2) = 11.2 km s 1 and 
b(H) = 15.9 km s" 1 . 
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Table 8. As Tableg] for fe(CO) = 0.3 km s"\ r ex (CO) = 5 K and Ar( 12 CO)//V( 13 CO) = 35. a 
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These shielding functions were computed for the same parameters as in Table|5] except A'( 12 CO)/A'( 13 CO) = 35. 
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Table 10. Column densities (in cm 2 ) for our grid of translucent cloud models at r gas = 15 K. Results for n H = 500 and 700 cm -3 
are not shown. 
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1.6(8) 


6.3(8) 


5.0(9) 


2.5(6) 


9.1(6) 








0.30 


1.0(20) 


2.8(20) 


2.5(12) 


1.2(9) 


4.6(9) 


3.7(10) 


1.8(7) 


6.8(7) 








0.70 


3.0(20) 


5.1(20) 


3.7(13) 


1.5(10) 


5.4(10) 


4.5(11) 


2.2(8) 


8.0(8) 








1.7 


1.0(21) 


6.7(20) 


4.4(14) 


1.3(11) 


4.5(11) 


1.9(13) 


5.8(9) 


2.0(10) 








4.2 


3.0(21) 


6.7(20) 


4.0(16) 


1.6(13) 


5.6(13) 


6.8(14) 


2.6(11) 


9.2(11) 
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Table 11. As Table [TU1 for T gas = 30 K. 



Model 




X 


A v 


2V(H 2 ) 


N(B) 


N( 12 


CO) 


N(C 17 0) 


iV(C lg O) 


N( 13 


CO) 


N( n C n O) 


Af( 13 C ls O) 




(cm ) 




(mag) 






















Bla 


100 


1 


0.06 


1.0(19) 


8.2(19) 


1.4( 


11) 


6.9(7) 


2.7(8) 


2.1 


9) 


1.0(6) 


3.8(6) 








0.12 


3.0(19) 


1.2(20) 


6.2 


11) 


3.0(8) 


1.2(9) 


9.1 


9) 


4.5(6) 


1.7(7) 








0.26 


1.0(20) 


2.2(20) 


4.1 


12) 


1.9(9) 


7.3(9) 


5.8 


10) 


2.9(7) 


1.1(8) 








0.64 


3.0(20) 


4.2(20) 


5.4 


13) 


2.0(10) 


7.4(10) 


6.0 


11) 


2.9(8) 


1.1(9) 








1.6 


1.0(21) 


5.8(20) 


6.2 


1 A \ 

14) 


1.7(11) 


5.9(11) 


1.1 


13) 


3.7(9) 


1.3(10) 








4 2 






D.'-r 


io; 


1 fin w 




1 ^ 

1 ,D 


1 SI 

i j) 


fi inn 


2.2(12) 


mp 


100 


3 


0.15 


1.0(19) 


2.2(20) 


3.2( 


10) 


1.6(7) 


5.9(7) 


4.6 


8) 


2.3(5) 


8.4(5) 








0.23 


3.0(19) 


3.0(20) 


1.5 


11) 


7.4(7) 


2.8(8) 


2.2 


9) 


1.1(6) 


4.1(6) 








0.40 


1.0(20) 


4.3(20) 


1.2 


12) 


5.8(8) 


2.2(9) 


1.8 


10) 


8.8(6) 


3.2(7) 








0.80 


3.0(20) 


6.6(20) 


1.7 


13) 


7.3(9) 


2.7(10) 


2.2 


11) 


1.1(8) 


4.0(8) 








1.8 


1.0(21) 


7.8(20) 


1.0 


14) 


3.7(10) 


1.3(11) 


2.0 


12) 


8.1(8) 


2.9(9) 








A Q 

4.3 


3.U(Z1 ) 


l.y(ZU) 


Z.t) 


16) 


y.4(lZ) 


1 Aft 1\ 

3.4(13) 


1.6 


1 A \ 

14) 


3.0(1 1 ) 


1 A fi 0\ 

1.4(12) 


Bly 


100 


10 


0.37 


1.0(19) 


5.6(20) 


1.1 


10) 


5.3(6) 


2.0(7) 


1.6 


8) 


7.9(4) 


2.9(5) 








0.48 


3.0(19) 


7.0(20) 


5.9 


10) 


2.8(7) 


1.1(8) 


8.7 


8) 


4.3(5) 


1.6(6) 








0.70 


1.0(20) 


9.0(20) 


5.8 


11) 


2.8(8) 


1.1(9) 


8.7 


9) 


4.3(6) 


1.6(7) 








1.1 


3.0(20) 


1.2(21) 


9.4 


12) 


4.3(9) 


1.6(10) 


1.3 


11) 


6.5(7) 


2.4(8) 








2.2 


1.0(21) 


1.4(21) 


1.5 


14) 


5.4(10) 


2.0(H) 


2.5 


12) 


1.0(9) 


3.7(9) 








A 1 

4. / 


3.0(Z1) 


1.4(21) 


1 r 
1 .J 


16) 


b.y( VI) 


L. j(1 j) 


J.y 


1 A \ 

14) 


Z.8(l 1 ) 


1 fV 1 o\ 


B2a 


300 


1 


0.03 


1.0(19) 


3.2(19) 


6.4( 


11) 


3.0(8) 


1.2(9) 


9.3 


9) 


4.6(6) 


1.7(7) 








0.07 


3.0(19) 


5.7(19) 


2.7 


12) 


1.3(9) 


4.9(9) 


3.9 


10) 


1.9(7) 


7.1(7) 








0.21 


1.0(20) 


1.3(20) 


1.7 


13) 


7.2(9) 


2.7(10) 


2.2 


11) 


1.1(8) 


4.0(8) 








0.56 


3.0(20) 


2.9(20) 


2.3 


14) 


6.7(10) 


2.4(11) 


2.1 


12) 


9.3(8) 


3.4(9) 








1 .5 


1.0(21) 


3.8(20) 


2.1 


15) 


4.8(1 1) 


1.5(12) 


5.3 


13) 


1.3(10) 


3.9(10) 








4. 




j.y\^-\j } 




1 1 ) 


1 4/141 


^ 1 f 1 4"! 
j. il.it; 


fi 


1 Si 

i j) 




s fin 7~\ 


B2/3 


300 


3 


0.06 


1.0(19) 


8.2(19) 


1.1 


11) 


5.4(7) 


2.1(8) 


1.6 


9) 


8.1(5) 


3.0(6) 








0.12 


3.0(19) 


1.2(20) 


5.4 


11) 


2.6(8) 


1.0(9) 


7.9 


9) 


3.9(6) 


1.4(7) 








0.26 


1.0(20) 


2.2(20) 


3.7( 


12) 


1.7(9) 


6.7(9) 


5.3 


10) 


2.6(7) 


9.8(7) 








0.64 


3.0(20) 


4.1(20) 


4.9 


13) 


1.9(10) 


6.8(10) 


5.6 


11) 


2.7(8) 


1.0(9) 








1.6 


1.0(21) 


5.8(20) 


6.K 


14) 


1.7(11) 


5.6(11) 


1.2 


13) 


3.6(9) 


1.2(10) 








4.2 


3.0(21) 


5.8(20) 


2.3 


17) 


9.2(13) 


1 A f \ A\ 

3.4(14) 


4.0 


15) 


l.o(lz) 


6.6(12) 


B2y 


300 


10 


0.16 


1.0(19) 


2.4(20) 


2.2( 


10) 


1.1(7) 


4.0(7) 


3.2 


8) 


1.6(5) 


5.8(5) 








0.24 


3.0(19) 


3.2(20) 


1.2 


11) 


5.6(7) 


2.2(8) 


1.7 


9) 


8.5(5) 


3.1(6) 








0.42 


1.0(20) 


4.7(20) 


1.0 


12) 


4.9(8) 


1.9(9) 


1.5 


10) 


7.5(6) 


2.8(7) 








0.82 


3.0(20) 


7.0(20) 


1.5 


13) 


6.5(9) 


2.4(10) 


2.0 


11) 


9.7(7) 


3.5(8) 








1.9 


1.0(21) 


9.2(20) 


2.3( 


14) 


7.4(10) 


2.6(H) 


3.9 


12) 


1.5(9) 


5.2(9) 








1 1 

4.4 


o.U(zl) 


y.Z(ZU) 


1 r 

1. J 


1') 


o,U(13) 


O 1/1 A\ 

2.2(14) 


1 A 

3.4 


15) 


1 K 1 1\ 

1.3(12) 


J.U(lZ) 


B5a 


1000 


1 


0.02 


1.0(19) 


1.4(19) 


3.6( 


12) 


1.7(9) 


6.5(9) 


5.1 


10) 


2.5(7) 


9.3(7) 








0.06 


3.0(19) 


3.1(19) 


1.5 


13) 


6.3(9) 


2.4(10) 


1.9 


11) 


9.5(7) 


3.5(8) 








0.18 


1.0(20) 


8.3(19) 


1.0 


14) 


3.2(10) 


1.2(11) 


1.0 


12) 


4.7(8) 


1.7(9) 








0.50 


3.0(20) 


1.9(20) 


1.3 


15) 


2.8(11) 


8.9(11) 


9.4 


12) 


3.5(9) 


1.3(10) 








1 .4 


1.0(21) 


2.0(20) 


3.2 


16) 


1 A/ 1 T\ 

o.U(12) 


o o/ i o\ 
o.o(lz) 


o. / 


1 A \ 

14) 


6.1(10) 


1.9(11) 








3.9 


3.0(21) 




5.5 


17^ 






8.5 


1 51 

i J ) 


2.6(12) 


9 7(1 7\ 
y. 


B5/J 


1000 


3 


0.03 


1.0(19) 


3.0(19) 


6.9 


11) 


3.3(8) 


1.3(9) 


1.0 


10) 


5.0(6) 


1.9(7) 








0.07 


3.0(19) 


5.4(19) 


3.0 


12) 


1.4(9) 


5.4(9) 


4.3 


10) 


2.1(7) 


7.9(7) 








0.20 


1.0(20) 


1.2(20) 


1.9( 


13) 


7.9(9) 


3.0(10) 


2.4 


11) 


1.2(8) 


4.4(8) 








0.56 


3.0(20) 


2.8(20) 


2.5 


14) 


7.3(10) 


2.6(11) 


2.3 


12) 


1.0(9) 


3.7(9) 








1.5 


1.0(21) 


3.7(20) 


2.8 


15) 


5.8(11) 


1.7(12) 


7.5 


13) 


1.6(10) 


4.8(10) 








4.0 


3.0(21) 


3.7(20) 


4.3 


17) 


1.4(14) 


5.1(14) 


6.8 


15) 


2.1(12) 


8.0(12) 


B5y 


1000 


10 


0.06 


1.0(19) 


8.2(19) 


1.0 


11) 


4.8(7) 


1.9(8) 


1.5 


9) 


7.3(5) 


2.7(6) 








0.12 


3.0(19) 


1.2(20) 


5.K 


11) 


2.4(8) 


9.4(8) 


7.5 


9) 


3.7(6) 


1.4(7) 








0.26 


1.0(20) 


2.2(20) 


3.6 


12) 


1.7(9) 


6.4(9) 


5.2 


10) 


2.6(7) 


9.5(7) 








0.64 


3.0(20) 


4.1(20) 


4.7( 


13) 


1.8(10) 


6.6(10) 


5.5 


11) 


2.6(8) 


9.7(8) 








1.6 


1.0(21) 


5.8(20) 


6.0 


14) 


1.7(11) 


5.5(11) 


1.2 


13) 


3.6(9) 


1.2(10) 








4.2 


3.0(21) 


5.8(20) 


3.4 


17) 


1.0(14) 


3.9(14) 


5.5 


15) 


1.6(12) 


6.2(12) 
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Table 12. As Table [TO] for 7/ gas = 50 K. 



Model 




X 


A v 


2V(H 2 ) 


N(B) 


N( 12 


CO) 


N(C 17 0) 


iV(C lg O) 


N( 13 


CO) 


N( n C n O) 


JV( 13 C ls O) 




(cm ) 




(mag) 






















Cla 


100 


1 


0.05 


1.0(19) 


6.5(19) 


1.9( 


11) 


9.3(7) 


3.6(8) 


2.8 


9) 


1.4(6) 


5.1(6) 








0.10 


3.0(19) 


1.0(20) 


8.5 


11) 


4.1(8) 


1-6(9) 


1.2 


10) 


6.2(6) 


2.3(7) 








0.24 


1.0(20) 


1.8(20) 


5.5 


12) 


2.5(9) 


9.6(9) 


7.7 


10) 


3.8(7) 


1.4(8) 








0.61 


3.0(20) 


3.6(20) 


6.9 


13) 


2.5(10) 


9.1(10) 


7.5 


11) 


3.6(8) 


1.3(9) 








1.5 


1.0(21) 


4.5(20) 


4.5 


14) 


1.3(11) 


4.4(1 1) 


7.6 


12) 


2.6(9) 


8.8(9) 








4 1 




4 fir?01 


1 2 


1 1 ) 






7 S 


1 J) 


1.2(12) 


4.2(12) 


c\p 


100 


3 


0.12 


1.0(19) 


1.7(20) 


4.0( 


10) 


1.9(7) 


7.4(7) 


5.8 


8) 


2.9(5) 


1.0(6) 








0.19 


3.0(19) 


2.4(20) 


2.0 


11) 


9.5(7) 


3.6(8) 


2.9 


9) 


1-4(6) 


5.2(6) 








0.36 


1.0(20) 


3.6(20) 


1.5 


12) 


7.3(8) 


2.8(9) 


2.2 


10) 


1.1(7) 


4.1(7) 








0.74 


3.0(20) 


5.8(20) 


2.0 


13) 


8.6(9) 


3.2(10) 


2.6 


11) 


1.3(8) 


4.7(8) 








1.7 


1.0(21) 


7.1(20) 


1.6 


14) 


5.5(10) 


1.9(11) 


2.5 


12) 


9.8(8) 


3.5(9) 








4.3 


J.U(21 ) 


1 .2(2U) 


H 1 


16) 


3.2(13) 


1 If i A\ 

1.2(14) 


1 .0 


15) 


/.5(ii) 


2. /(12) 


Cly 


100 


10 


0.31 


1.0(19) 


4.6(20) 


1.2 


10) 


5.7(6) 


2.2(7) 


1.7 


8) 


8.4(4) 


3.1(5) 








0.41 


3.0(19) 


6.0(20) 


6.5 


10) 


3.1(7) 


1.2(8) 


9.6 


8) 


4.7(5) 


1.7(6) 








0.62 


1.0(20) 


7.8(20) 


6.4 


11) 


3.1(8) 


1.2(9) 


9.5 


9) 


4.7(6) 


1.7(7) 








1.0 


3.0(20) 


1.0(21) 


9.8 


12) 


4.5(9) 


1.7(10) 


1.4 


11) 


6.8(7) 


2.5(8) 








2.0 


1.0(21) 


1.2(21) 


5.5( 


13) 


2.1(10) 


7.9(10) 


8.5 


11) 


3.7(8) 


1.3(9) 








4. 3 


7 n/"o 1 \ 
J.U(Zl) 


1.2(21) 


4. j 


16) 


2.1(13) 


/.4(13) 


t i 


15) 


5.1(11) 


1 Qi 1 T\ 

1.0(12) 


C2a 


300 


1 


0.03 


1.0(19) 


2.6(19) 


9.K 


11) 


4.3(8) 


1-7(9) 


1.3 


10) 


6.5(6) 


2.4(7) 








0.07 


3.0(19) 


4.8(19) 


3.8 


12) 


1.8(9) 


6.8(9) 


5.4 


10) 


2.7(7) 


9.9(7) 








0.20 


1.0(20) 


1.1(20) 


2.4 


13) 


9.6(9) 


3.6(10) 


2.9 


11) 


1.4(8) 


5.3(8) 








0.54 


3.0(20) 


2.6(20) 


3.0 


14) 


8.5(10) 


3.0(11) 


2.6 


12) 


1.2(9) 


4.3(9) 








1 .5 


1.0(21) 


3.4(20) 


3.2 


15) 


6.3(1 1) 


1.9(12) 


5.5 


13) 


1.2(10) 


3.7(10) 








4 






j.y 


1 1 ) 






u.z. 


1 

i j) 


? 3f19"l 




C2J3 


300 


3 


0.05 


1.0(19) 


6.5(19) 


1.6 


11) 


7.8(7) 


3.0(8) 


2.4 


9) 


1.2(6) 


4.3(6) 








0.10 


3.0(19) 


1.0(20) 


7.7 


11) 


3.7(8) 


1.4(9) 


1.1 


10) 


5.6(6) 


2.1(7) 








0.24 


1.0(20) 


1.8(20) 


5.1 


12) 


2.4(9) 


9.0(9) 


7.2 


10) 


3.6(7) 


1.3(8) 








0.61 


3.0(20) 


3.6(20) 


6.5 


13) 


2.3(10) 


8.6(10) 


7.1 


11) 


3.4(8) 


1.3(9) 








1.5 


1.0(21) 


4.5(20) 


3.8 


14) 


1.1(11) 


3.8(11) 


6.2 


12) 


2.1(9) 


7.3(9) 








i i 
4. 1 


3.0(21) 


4.5(20) 


2.8 


17) 


1 1 (\ A \ 

1.1(14) 


1 C\f 1 A \ 

3.^(14) 


4.0 


15) 


1 Hi 1 0\ 

1.7(12) 


0.4(12) 


C2y 


300 


10 


0.13 


1.0(19) 


1.9(20) 


2.9( 


10) 


1.4(7) 


5.3(7) 


4.2 


8) 


2.1(5) 


7.6(5) 








0.20 


3.0(19) 


2.6(20) 


1.6 


11) 


7.5(7) 


2.9(8) 


2.3 


9) 


1.1(6) 


4.2(6) 








0.37 


1.0(20) 


3.9(20) 


1.3 


12) 


6.3(8) 


2.4(9) 


1.9 


10) 


9.6(6) 


3.5(7) 








0.76 


3.0(20) 


6.1(20) 


1.8 


13) 


7.7(9) 


2.9(10) 


2.3 


11) 


1.1(8) 


4.2(8) 








1.7 


1.0(21) 


7.3(20) 


1.1 


14) 


4.0(10) 


1.5(11) 


1.7 


12) 


7.1(8) 


2.6(9) 








4. J 


3.U(21) 


i.i(ZU) 


1 o 


l/) 


l.vyli) 


2.0(14) 


i.i 


15) 


1.2(12) 


4.0(12) 


C5a 


1000 


1 


0.02 


1.0(19) 


1.2(19) 


5.0( 


12) 


2.2(9) 


8.7(9) 


6.8 


10) 


3.4(7) 


1.2(8) 








0.05 


3.0(19) 


2.7(19) 


2.K 


13) 


8.3(9) 


3.2(10) 


2.6 


11) 


1.3(8) 


4.6(8) 








0.17 


1.0(20) 


7.2(19) 


1.4 


14) 


4.2(10) 


1.6(11) 


1.3 


12) 


6.1(8) 


2.2(9) 








0.48 


3.0(20) 


1.7(20) 


1.8 


15) 


3.7(11) 


1.1(12) 


1.2 


13) 


4.4(9) 


1.6(10) 








1 .4 


1.0(21) 


2.0(20) 


O 1 

L.i 


16) 


2.3(12) 


6.6(12) 


4.2 


1 A \ 

14) 


4.2(10) 


1 Iff 1\ 

1.3(11) 








3.9 


3.0(21) 




5.4 


17^ 




6.2(14) 


8.0 


1 

1 J ) 


2.5(12) 


9 3(1 7 s ! 


C5/3 


1000 


3 


0.03 


1.0(19) 


2.4(19) 


1.0 


12) 


4.7(8) 


1.8(9) 


1.5 


10) 


7.2(6) 


2.7(7) 








0.07 


3.0(19) 


4.5(19) 


4.2 


12) 


1.9(9) 


7.6(9) 


6.0 


10) 


3.0(7) 


1.1(8) 








0.19 


1.0(20) 


1.1(20) 


2.7 


13) 


1.1(10) 


4.0(10) 


3.3 


11) 


1.6(8) 


5.9(8) 








0.54 


3.0(20) 


2.5(20) 


3.4 


14) 


9.4(10) 


3.2(11) 


2.9 


12) 


1.3(9) 


4.7(9) 








1.5 


1.0(21) 


3.3(20) 


4.0 


15) 


7.2(11) 


2.1(12) 


7.1 


13) 


1.4(10) 


4.2(10) 








4.0 


3.0(21) 


3.4(20) 


4.4 


17) 


1.4(14) 


5.1(14) 


6.7 


15) 


2.0(12) 


7.7(12) 


C5y 


1000 


10 


0.05 


1.0(19) 


6.5(19) 


1.5 


11) 


7.3(7) 


2.8(8) 


2.2 


9) 


1.1(6) 


4.0(6) 








0.10 


3.0(19) 


1.0(20) 


7.4( 


11) 


3.5(8) 


1.4(9) 


1.1 


10) 


5.4(6) 


2.0(7) 








0.24 


1.0(20) 


1.8(20) 


5.0 


12) 


2.3(9) 


8.8(9) 


7.1 


10) 


3.5(7) 


1.3(8) 








0.61 


3.0(20) 


3.6(20) 


6.4 


13) 


2.3(10) 


8.4(10) 


7.0 


11) 


3.4(8) 


1.2(9) 








1.5 


1.0(21) 


4.5(20) 


3.5 


14) 


1.0(11) 


3.5(11) 


5.7 


12) 


1.9(9) 


6.8(9) 








4.1 


3.0(21) 


4.5(20) 


3.3 


17) 


9.8(13) 


3.7(14) 


5.1 


15) 


1.5(12) 


5.7(12) 
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Table 13. As Table [TO] for 7/ gas = 100 K. 



Model 




X 


A v 


2V(H 2 ) 


2V(H) 


N( 12 


CO) 


N(C 17 0) 


iV(C lg O) 


N( 13 


CO) 


N( n C n O) 


Af( 13 C ls O) 




(cm J 




(mag) 






















Dior 


100 


1 


0.04 


1.0(19) 


4.8(19) 


3.5( 


11) 


1.6(8) 


6.4(8) 


5.0 


9) 


2.5(6) 


9.1(6) 








0.09 


3.0(19) 


7.8(19) 


1.5 


12) 


7.1(8) 


2.8(9) 


2.2 


10) 


1.1(7) 


4.0(7) 








0.22 


1.0(20) 


1.6(20) 


9.5 


12) 


4.2(9) 


1.6(10) 


1.3 


11) 


6.3(7) 


2.3(8) 








0.59 


3.0(20) 


3.3(20) 


1.2 


14) 


3.9(10) 


1.4(11) 


1.2 


12) 


5.5(8) 


2.0(9) 








1.5 


1.0(21) 


4.2(20) 


8.7 


14) 


2.3(11) 


7.4(11) 


1.2 


13) 


3.6(9) 


1.2(10) 








4 




A onm 

^t.Z^ZU f 


2.7 


1 71 






4 1 






fi fin ?i 


DlyS 


100 


3 


0.09 


1.0(19) 


1.3(20) 


6.7 


10) 


3.2(7) 


1.2(8) 


9.6 


8) 


4.8(5) 


1.8(6) 








0.15 


3.0(19) 


1.8(20) 


3.4 


11) 


1.6(8) 


6.3(8) 


5.0 


9) 


2.4(6) 


9.0(6) 








0.31 


1.0(20) 


2.9(20) 


2.5 


12) 


1.2(9) 


4.5(9) 


3.6 


10) 


1.8(7) 


6.6(7) 








0.69 


3.0(20) 


5.0(20) 


3.2 


13) 


1.3(10) 


4.8(10) 


3.9 


11) 


1-9(8) 


6.9(8) 








1.7 


1.0(21) 


6.5(20) 


2.9 


14) 


9.0(10) 


3.1(11) 


3.5 


12) 


1.4(9) 


4.9(9) 








A 1 

4.z 


i.U(Zl ) 


O.j(zU) 


1 a 

1 .0 


1') 


/.U(13) 


2.5(14) 


Z. 1 


15) 


l.z(lZ) 


4.2(12) 


Dly 


100 


10 


0.24 


1.0(19) 


3.6(20) 


1.6 


10) 


7.8(6) 


3.0(7) 


2.3 


8) 


1.1(5) 


4.2(5) 








0.33 


3.0(19) 


4.7(20) 


9.4 


10) 


4.5(7) 


1.7(8) 


1.4 


9) 


6.8(5) 


2.5(6) 








0.53 


1.0(20) 


6.3(20) 


9.1 


11) 


4.4(8) 


1.7(9) 


1.3 


10) 


6.7(6) 


2.4(7) 








0.93 


3.0(20) 


8.7(20) 


1.3 


13) 


5.9(9) 


2.2(10) 


1.8 


11) 


8.9(7) 


3.2(8) 








1.9 


1.0(21) 


1.0(21) 


1.1 


14) 


3.9(10) 


1.4(11) 


1.4 


12) 


6.0(8) 


2.2(9) 








1 c 

4. J 


J.O(zl) 


1 1 /"") 1 \ 

1.1(21) 


1 .U 


i ') 


4.j(13) 


1.0(14) 


1.5 


15) 


/./(ll) 


z.o(lz) 


D2a 


300 


1 


0.03 


1.0(19) 


2.0(19) 


1.6 


12) 


7.6(8) 


3.0(9) 


2.3 


10) 


1.2(7) 


4.3(7) 








0.06 


3.0(19) 


4.0(19) 


6.7 


12) 


3.0(9) 


1.2(10) 


9.2 


10) 


4.6(7) 


1.7(8) 








0.19 


1.0(20) 


9.8(19) 


4.3 


13) 


1.6(10) 


5.9(10) 


4.8 


11) 


2.3(8) 


8.6(8) 








0.52 


3.0(20) 


2.3(20) 


5.2 


14) 


1.3(11) 


4.5(11) 


4.1 


12) 


1.8(9) 


6.4(9) 








1 .5 


1.0(21) 


3.1(20) 


6. 1 


15) 


9.6(1 1) 


2.8(12) 


7.8 


13) 


1.4(10) 


4.5(10) 








4 


^ (¥911 


j.z^zu j 


4. 1 


i 1 ) 


1 fin 4.1 




fi A 


1 •>) 


? 3fiyi 




D2/J 


300 


3 


0.04 


1.0(19) 


4.8(19) 


3.1 


11) 


1.5(8) 


5.8(8) 


4.6 


9) 


2.3(6) 


8.3(6) 








0.09 


3.0(19) 


7.8(19) 


1.4 


12) 


6.8(8) 


2.6(9) 


2.1 


10) 


1.0(7) 


3.8(7) 








0.22 


1.0(20) 


1.5(20) 


9.1 


12) 


4.1(9) 


1.5(10) 


1.2 


11) 


6.1(7) 


2.3(8) 








0.59 


3.0(20) 


3.2(20) 


1.1 


14) 


3.8(10) 


1.4(11) 


1.1 


12) 


5.4(8) 


2.0(9) 








1.5 


1.0(21) 


4.2(20) 


6.6 


14) 


1.7(11) 


5.8(11) 


8.1 


12) 


2.7(9) 


9.3(9) 








4.0 


3.0(21) 


4.2(20) 


1 o 

i.l 


1 ~7\ 

17) 


1.2(14) 


4.3(14) 


4.8 


15) 


l.o(lz) 


0.4(12) 


D2y 


300 


10 


0.10 


1.0(19) 


1.4(20) 


5.2 


10) 


2.5(7) 


9.5(7) 


7.5 


8) 


3.7(5) 


1.4(6) 








0.16 


3.0(19) 


2.0(20) 


2.8 


11) 


1.4(8) 


5.2(8) 


4.1 


9) 


2.0(6) 


7.5(6) 








0.32 


1.0(20) 


3.1(20) 


2.2 


12) 


1.0(9) 


4.0(9) 


3.2 


10) 


1.6(7) 


5.8(7) 








0.71 


3.0(20) 


5.2(20) 


2.8 


13) 


1.2(10) 


4.3(10) 


3.5 


11) 


1.7(8) 


6.2(8) 








1.7 


1.0(21) 


6.6(20) 


2.0 


14) 


6.4(10) 


2.3(H) 


2.3 


12) 


9.7(8) 


3.5(9) 








4.z 


J.U(zl) 


0. /(ZU) 


z.z 


1 ') 


Q 1 /1 1\ 
O. 1(1 J) 


o.0(14) 




15) 


1.2(12) 


A A/ 1 T\ 
4.0(lz) 


D5o- 


1000 


1 


0.02 


1.0(19) 


9.7(18) 


8.4 


12) 


3.6(9) 


1.4(10) 


1.1 


11) 


5.4(7) 


2.0(8) 








0.05 


3.0(19) 


2.3(19) 


3.6 


13) 


1.3(10) 


5.1(10) 


4.1 


11) 


2.0(8) 


7.3(8) 








0.17 


1.0(20) 


6.3(19) 


2.4 


14) 


6.6(10) 


2.4(11) 


2.1 


12) 


9.2(8) 


3.4(9) 








0.47 


3.0(20) 


1.4(20) 


3.2 


15) 


5.8(11) 


1.7(12) 


1.9 


13) 


6.6(9) 


2.3(10) 








1 .4 


1.0(21 ) 


1. /(20) 


Q a 

3.6 


16) 


J.U(lz) 


0.5(12) 


A "7 
4. / 


1 1 \ 

14) 


4.4(10) 


1.4(11) 








3.9 


3.0(21) 


1.7(20) 


5.6 


171 




6.2(14) 


8.1 


1 51 


2.4(12) 


9 1C191 


D5/3 


1000 


3 


0.02 


1.0(19) 


1.9(19) 


1.8 


12) 


8.6(8) 


3.4(9) 


2.6 


10) 


1.3(7) 


4.8(7) 








0.06 


3.0(19) 


3.8(19) 


7.6 


12) 


3.4(9) 


1.3(10) 


1.0 


11) 


5.1(7) 


1.9(8) 








0.19 


1.0(20) 


9.4(19) 


4.9 


13) 


1.8(10) 


6.6(10) 


5.4 


11) 


2.6(8) 


9.6(8) 








0.52 


3.0(20) 


2.2(20) 


5.9 


14) 


1.5(11) 


5.0(11) 


4.6 


12) 


2.0(9) 


7.1(9) 








1.5 


1.0(21) 


3.0(20) 


7.3( 


15) 


1.0(12) 


3.0(12) 


9.1 


13) 


1.5(10) 


4.8(10) 








4.0 


3.0(21) 


3.0(20) 


4.6( 


17) 


1.4(14) 


5.1(14) 


6.8 


15) 


2.0(12) 


7.5(12) 


D5y 


1000 


10 


0.04 


1.0(19) 


4.8(19) 


3.0 


11) 


1.5(8) 


5.6(8) 


4.4 


9) 


2.2(6) 


8.1(6) 








0.09 


3.0(19) 


7.8(19) 


1.4 


12) 


6.7(8) 


2.6(9) 


2.0 


10) 


1.0(7) 


3.7(7) 








0.22 


1.0(20) 


1.5(20) 


9.0 


12) 


4.0(9) 


1.5(10) 


1.2 


11) 


6.0(7) 


2.2(8) 








0.58 


3.0(20) 


3.2(20) 


1.1 


14) 


3.7(10) 


1.3(11) 


1.1 


12) 


5.3(8) 


1.9(9) 








1.5 


1.0(21) 


4.2(20) 


5.8 


14) 


1.6(11) 


5.2(11) 


6.9 


12) 


2.4(9) 


8.3(9) 








4.0 


3.0(21) 


4.2(20) 


3.4 


17) 


9.8(13) 


3.7(14) 


5.0 


15) 


1.4(12) 


5.5(12) 



